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FOREWORD 


This document describes the work performed during Phase 
VII of JPL Contract 952532. (This is a follow-on to JPL 
Contract 95202-8). Phase VII accomplished improvements in the 
computer programs generated in Phase III. 

This document also contains a revised User's Manual for 
the Microbial Burden Prediction Model and supplants the pre- 
vious User's Manual (Volume II) for this purpose. Volume II 
is still useful for details of the theory of the model. 
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DEFINITIONS 


Alphameric - (Short for Alphabetic-numeric) 

Data such as "PART 2W31" that may consist of a 
mixture of letters, numerals, and punctuation. 

Card Image The data on a tape or disk that represents a 
punched card. 


Characte r The contents of one column on a punched card or 
its equivalent on tape or disk. 


Disk A storage device in a computer. 

Field A group of columns on a punched card that is devoted 

to one quantity. 


Float inp: Point 

A number expressed as a characteristic and an exponent 

4 

as in scientific notation: 3.2 x 10 . When punched 

on cards or printed by the computer, the above number 
would appear as 3.2E+4. (See Integer. ) 

Histogram (Used interchangeably with "distribution" in this report). 

A method of representing a probability distribution by 
means of a number of intervals with a probability asso- 
ciated with each; a "bar graph" of the probability density. 

Integer A computer quantity opposed to Floating Point. An integer 

or "fixed point" number repsents a whole number, has no 
decimal point, and must be as far right in its field as 
it can go (right justified) on a punched card. 


Record See Tape 

Stage, Task, Subtask, Operation - Terms used to denote the level of 

activity in an assembly and test sequence to be simulated 
by the burden prediction model. The operation is the most 
detailed level. 
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Tape 


Word 


Definitions (Continued) 

Magnetic tape, a storage medium used by computers. 
Information on tapes (or disks) consists of "logical 
records" each of which contains a number of computer 
words. Information is written on and read off tapes 
(or disks) by logical records. 

A storage location in a computer that can contain either 
a single number (integer or floating point) or several 
alphameric characters. 


v 



I. 


INTRODUCTION 


During Phases I, II, and III of JPL Contract 952028, the requirements 
for a microbial burden prediction model were established, and the model was 
developed and programmed for computer., During Phases IV, V, and VI, the 
model was used to estimate the burden on the Mariner V Spacecraft* Refer- 
ences 1 through 5 describe this work. 

As a result of the above work, a number of changes were identified that 

would improve the effectiveness of the existing computer programs. These 

changes were grouped into four tasks (JPL Contract 952532) as described in 
the following sections. 

A. Daily Restarts . This modification allows new data (e.g. , the assembly 
and test activities performed during one day) to be incorporated to update a 
previous estimate. The results of the previous run are written on tape and 
saved to initialize the new run. This capability permits frequent updates 
without the wasted computer time involved in recomputing from the entire data 
deck. An advantage is that, if desired, a tape may be Whitten for the last 
task successfully completed so that an abnormal job termination v/ill not 
require recomputing the entire run. 

In making this change, the original restart subroutine was removed and 
all "write” statements for Tape 11 (which prevously had recorded both restart 
and detailed data) were modified to contain only the detailed information not 
normally printed during a run. A data retrieval program was prepared to print 
information on selected tasks and subtasks from Tape 11. 



B. Alteration of Data on Tape . The input data cards for the prediction 
program were previously read and checked by a data check program and then 
written on a tape (Tape 9)* Any subsequent changes in data required changing 
the data cards and re-running the data check program. To avoid working with 
the card deck, a computer program was prepared to alter the data already on 
tape, including checking the resulting data and writing it on a new tape. 
Changes possible by this method consist of: 

a) Changing any parameter on any card; 

b) Adding new data (from cards) at any point in the tape; 

c) Changing any parameter (either by replacing it or by 
multiplying it by a constant) every time it occurs 
between specified points on the tape; 

d) Deleting specified card images from the tape. 

The capability b) also permits data inputs entirely from cards in case no 
data tape is used. 

As a data tape is prepared, a listing is made identifying each tape record 
by reference number (sequential), record type, stage, task, and subtask. The 
reference number is used for specifying future changes in that tape. 

C. Faster Running . Burden calculations are performed using histograms, and 
the time required to perform a given run varies approximately as the square 

of the number of histogram intervals used. The capability was added to permit 
specifying the number of intervals (up to ten) to be used for each task. This 
allows a rough estimation to be run with one or two intervals (or even no 
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intervals; that is using mean values) using little computer time. 

D. Print-out Options . The Mariner V runs took over two hours to print out 
the results. To reduce paper volume and print-out time, a more compact out- 
put format is now used. This output can be specified to be at the subtask, 
operation, or part level as desired for each task. Burden histograms are no 
longer printed except for inputs and task summaries (all burden histograms 
are recorded on Tape 11 for printing if desired) , but some information has 
been included in the normal printout that was not printed before (e.g., the 
number of men working near the hardware). 


During these four tasks, it was found useful to combine the programs for 
tape alteration (including card reading), burden estimation, and detailed 
printout from Tape 11 into one program so that separate submittals were not 
necessary. The combination requires more central memory storage in the com- 
puter than the individual programs do (unless segmentation or overlay is used). 
The three programs may easily be separated if the combination is not desired. 

The programs are designed to run on the JPL Univac 1108 computer, but 
were programmed so that they could run on other large computers with essentially 
no changes. 

The present volume (Volume VI) is complete as far as users of the computer 
programs are concerned; it supercedes all user instructions in the previous 
User’s Manual (Volume II). However, Volume II contains explanations of the 
model and the meanings of the variables used. Volume II should therefore be 
regarded as a reference rather than a user’s manual. 



II. TECHNICAL DISCUSSION (Revised User’s Manual) 


A. Descriptions of Changes Made 

1. Daily Restarts . To provide an improved restart capability, a 
dump is made (on Tape 2) of all information "remembered" in the computer 
from one control card to the next. This information is approximately that 
contained in blank common storage, but differs to the extent that a few 
extra variables are included and any unused histograms are omitted. The 
tape records are in binary, and will vary in number 100 to 1000 depending 
on the number of histograms used. 

Three options are provided for this feature: 

1) Restart dump at end of run only; 

2) Restart dump at the end of the last 
task successfully completed; 

3) No restart dumps. 

These options are selected by setting the variable 12 on any control card 
and can be changed during the run if desired. (See Section IIC for card 
preparation. ) Option #2 preserves computed results task by task so that 
an error exit from the program will require recomputation only from 
last completed task instead of re-running the entire data deck. 

The burden prediction program previously included a restart capability, 
embodied in subroutine MBRS, the Microbial Burden Restart Subroutine, which 
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made use of detailed information recorded on Taps. 11 during a previous run. 
Subroutine MBRS is no longer used, and Tape 11 is now used only for recording 
detailed burden information which is not normally of interest and would thus 
cause unnecessary printout volume (and time) if included in the printed output 
of the burden prediction program. 

All records written on Tape 11 are identical, consisting of 43 words. 

In general, each record contains one histogram, the run, stage, task, and sub- 
task numbers, the subtask description, and ten additional parameters. Each 
record is identified as to type, so that the information recorded in each word 
can be identified and printed. A separate computer program (subroutine DPS, 

PLS, and HPS) has been written to locate and print out selected tasks and sub- 
tasks from Tape' 11. The information to be printed is determined by a set of 
instruction cards, which can specify either a task (or subtask) to be printed 
or that all information between a given task (or subtask) and another. Any 
number of such instruction cards may be used. Specification of the subtask is 
optional; if not supplied, selection will be by tasks. See Table l'_ for typical 
output from this program. Table 2 describes the instruction types. 

To facilitate frequent updates in burden prediction, provision has been 
made for data inputs to the prediction program from cards instead of tape. 

The card data wse previously processed by a Data Check Program and written on 
tape (Tape 9 ) for later use by the prediction program, thus requiring two job 
submittals to obtain a burden prediction. The functions of data checking and 
recording (now accomplished by the Tape Alteration Subroutine; see Section II.A.2) 
have been combined in one program with the burden prediction program and the 



KA NGE 


1 . 58t-u 3 9.04E-t-02 


PART BURDEN AFTER UPDATE RUN 1, STAGE 7, TASK 5, SUBTASK 3 INST 800ST DAMPERS 

PART 5 SQX-3762-599, SURF. 1, AREA = .50, TIME= 30.000, ENVIRON. -0 BURDEN MEAN VALUE = 4.52129E+02 

PROBABILITY = 1.00000 

RANGE = 1.S8E-03 9.Q4E+02 


PART BURDEN AFTER DECONTAMINATION RUN 1, STAGE 7 , TASK 5, SUBTASK 8 INST BOOST DAMPERS 

PART 5, SURF. 1, AREA= .50, OPER. 1, .50 HOURS, FRACTION REMO VED= .500 BURDEN MEAN VALUE = 2.260&4E+02 

PROBABILITY = 1.00000 

RANGE = 3.94E-Q4 4.52E+02 


PART BURDEN AFTER FALLOUT RUN 1, STAGE 7 , TASK 5, SUSTASK 8 INST BOOST UAMPERS 

PART 5, SURF. 1, AREA = .50, OPER. 1, 10.00 HOURS, ENVIRON. 1, 1 MEN. BURDEN MEAN VALUE = 4.6G786E+G2 

PROBABILITY = 1.00000 

RANGE = 2.20E-04 9.22E+Q2 


PART BURDEN AFTER CONTACT RUN 1, STAGE 7 , TASK 5, SUBTASK 8 INST BOOST UAMPERS 

PART 5, SURF. 1, AREA= .50, CONTACT* .100, S1=7.0Q0, S2=1.QQ0, TOOL* 1.0 BURDEN MlAN VALUE = 4.15057E+02 

PROBABILITY = 1.00000 

RANGE = 2.1S3E-04 8.3QE+02 


START OF NEW SUBTASK RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PREREQUISITE SUBTASKS- -0, -0 START TIME MEAN VALUE = 2.64905t+C2 

PROBABILITY = 1.00000 

RANGE * 3.01E+Q1 5.00E+02 


PART BURDEN AFTER UPDATE RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PART 1 SQX-3755-420 , SURF. 1, ARcA* .10, TIME* 264.905, ENVIRON. -0 BURDEN MEAN VALUE = 3.69665E+02 

PROBABILITY = 1.00000 

RANGE = 6.09E-10 7.39E+02 


PART BURDEN AFTER FALLOUT RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PART 1, SURF. 1, AREA= .10, OPER. 1, 10.00 HOURS, ENVIRON. 1, 1 MEN. BURDEN MEAN VALUE = 3.66502E+02 

PROBABILITY * .10000 .20000 .30000 .15000 .10000 .05000 .04000 .03000 .02000 .01000 

RANGE = 3.58E-08 2.13E+01 6.40E+01 1.30E+02 2.16E+02 6.09E+02 8.05E+02 9.62E+02 1.08E+03 1.16E4Q3 2.62E+04 


PART BURDEN AFTER UPDATE RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PART 5 SQX-3762- 599 , SURF. 1, AREA* .50, T IME= 274.905, ENVIRON. -0 BURDEN MEAN VALUE = 4.15057E+02 

PROBABILITY * 1.00000 

RANGE = 2.10E-04 8.30E+02 


PART BURDEN AFTER FALLOUT RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PART 5 , SURF. 1, AREA= .50, OPER. 1, 10.00 HOURS, ENVIRON. 1, 3 MEN. BURDEN MEAN VALUE = 7.52761E+02 

PROBABILITY = .10000 .20000 .30000 .15000 .10000 .05000 .04000 .03000 .02000 .01000 

RANGE = 1.82E-07 9.66E+-00 2.75E+01 5.41E + 01 1.20E*02 1.97E+02 2.36E + 02 2.67E+02 2.9QE+02 8.98E*-Q2 1.32E + 05 

O' 

PART BURDEN AFTER CONTACT RUN 1, STAGE 7, TASK 99, SUBTASK 1 INSERT GUIDE PIN 

PART 5, SURF. 1, AREA* .50, CONTACT= .100, Sl=7.030, S2=1.000, TOOL= 1.0 BURDEN MEAN VALUE = 6.77335E<-02 

PROBABILITY = .10000 .20000 .30000 .15000 .10000 .05000 .04000 .03000 .02000 .01000 

RANGE = 7.&9E-06 5.47E+00 2.15E+01 4.64E+01 1.07E+02 1.81E+0Z 2.15E+02 2.51E+02 2.65E+02 5.2QE4Q2 1.20E+05 


Table 1. Typical Output, Detailed Printout Program 



Instruction Format: 


Columns: 5 9-10 13-15 19-20 24-25 28-30 34-35 

K NS NT NST NS2 NT2 NST2 

All numbers must be right justified in their fields. 

K determines the printout: 

K = 0: No more printing (program exit). 

K = 1: Print all information on stage NS, task NT, Subtask NST 

K = 2: Print all information from stage NS, task NT, subtask NST 

to (but not including) stage NS2, task NT2, subtask NST2. 

The subtask specification is optional; if omitted, the entire task 
will be printed. If all information beyond a specified task is 
desired, use a card with K = 2 and NS2 blank. 
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program to print from Tape 11 so that multiple job submittals are avoided. 

This was accomplished by making subroutines of the three programs and writing 
one main program that does nothing but call the three sub-programs in turn. 

If all three are not desired, the main program is easily altered (by removing 
appropriate "CALL..." cards) to accomplish this. The data interfaces between 
these sub-programs (Tape 9 between card reading and burden prediction, and 
Tape 11 between burden prediction and detailed printout) can be disc files 
unless permanent records are desired on tape; this option is determined by 
the presence or absence of tape request cards in the job control deck, the 
disc being used if a tape is not supplied. 

2. Alteration of Data on Tape . To alter a parameter on a card record on 
tape, three items of information are necessary: 

1) which tape record; 

2) which parameter on that record; 

3) the new value of the parameter. 

Tape records were constructed that, with one exception, are of uniform length: 
seven integer words plus seven floating point or alphameric words. The excep- 
tion is the histogram record, which consists of one integer plus 22 floating 
point words. The first word of each record identifies the record type and 
hence, the information contained in it. Table 3 shows the information layout 
of these records, and also shows the numbering system (from 1 through 13 except 
for histogram records, which are 1 through 22) used to identify the parameter 


to be changed. 



Table 3 Tape 9 Record Formats 


PARAMETERS 


Card 

Type 


J NTEGER 


FLOATING OR ALPHAMERIC 


CC 

RD 

SD 

TD 

EM 

ED 

EQ 

OD 

OQ 

PD 

PQ 

DD 

KD 


DQ 


3 


10 


_L2_ 


_L3 


KK 

KR 

KS 

KT 

N 

I EC 

N 

IOT 

N 


IK 


1 f 


12 


15 


1 6 


13 


14 




* 4 - 


L 1 


L2 


L3 


L4 


IET 


I OQ 


IOC 


|AB(1) I AB ( 2 ) I AB C 3 ) IABC4) !AB(5) !AB(6) 


M 

N I 


L 

N2 


KC 

K 

IRC!) 1 R (2 ) 

IRC3) 

IR(4) 

ARM ) 

AR(2) 

AR ( 3 ) 


KO 

10 

1 KE 



AKT 

AKQ 



PE 

1 P 

LS 



APD 

A PC 

APS 

ARAM ) 

PE 

1 P 

LK 



ARM ) 

ARC2) 

AR ( 3 ) 

AR(4) 

ZD 

IZ 




-< 


DSC (4) — 

5.. 

ZC 

1 P 




FPM ) 

FP(2) 

FP(3) 

FP(4) 


- RUN C 7 ) 

- STG (7) 
TSK (7) 








AES ( I ) AES (2) AES (3) AES(4) 


AEC 


AET 




DSCC4) 

aed 

DSC(4) 
DSC (6) 




AEF(1) AEF(2) AEF(3) AEFC4) 

a*. 


AAG 


AAS 




- DSC (4) 

- DSC(4) 




APAC2) 


Different Format 22 Parameters; DR(M>, DR(N,2) , . . . ,DR(N, 1 1 ) , XR(N,1). ,-XR(N 11) Numbered 1 thru 22. 


VO 





Identification of the tape record containing the desired parameter is 
most conveniently done by specifying the stage, task, subtask, etc. in the 
data deck since this is the way the data is organized. However duplicate 
task and subtask numbers are permitted, so such identification is not unique. 
For this and other reasons, the primary means of identifying the tape record 
to be altered is a single number, the number of that record in a sequential 
listing of all records on the tape. 

When a data tape (Tape 9 ) is written, a listing is made of the data on 
each record in the format of Table in addition to the record type and 
parameter values, the record number and the stage, task and subtask are 
printed. A typical listing is given in Table 4. 

The types of changes possible with this program are listed in Table 5« 
Essentially, any change that could be made in a card deck can be made in the 
tape records; the main difference is that it is not possible to alter just 
one character in tape record - the entire parameter must be replaced (or in 
the case of alphameric fields, a new card must be read to replace the record 
to be altered). 

When this program is used, the old data tape (if any) is designated 
Tape 12 and a new Tape 9 is written on a blank tape. A card deck must be 
provided to specify the changes to be made and to provide any new data cards 
to be added. Note that Tape 12 is not necessary; input can be entirely from 


cards if desired. 



RECORD, TYPE, S/TSK/ST 


2 3 4 5 6 7 


9 10 11 12 13 


CHANGE 


2 9999 


- 0 . 


1 

CC 



2 

RO 



3 

SO 

1/ 


4 

TO 

1/ 

6/ 

5 

EM 

1/ 

6/ 

6 

ED 

1/ 

6/ 

7 

EQ 

1/ 

6/ 

8 

ED 

1/ 

6/ 

9 

OO 

1/ 

6/ 

10 

OQ 

1/ 

6/ 

11 

OD 

1/ 

6/ 

12 

OQ 

1/ 

6/ 

13 

OD 

1/ 

6/ 

14 

PD 

1/ 

6/ 

15 

PQ 

1/ 

6/ 

16 

PD 

1/ 

6/ 

17 

DO 

1/ 

6/ 

18 

DQ 

1/ 

6/ 



O . 



-0. 


-0, 

19 

OD 

1/ 

6/ 

20 

UQ 

1/ 

6/ 



0. 



1.00E+00 

2 

21 

OD 

1/ 

6/ 

22 

DQ 

1/ 

6/ 



. 10 


‘'“'■ERROR 

1 

TOT/ 


l.OOE+OO 

3 

23 

DD 

1/ 

6/ 

24 

KD 

1/ 

6/ 

25 

KC 

1/ 

6/ 

26 

KO 

1/ 

6/ 

27 

PE 

1/ 

6/ 

28 

PE 

1/ 

6/ 

29 

KO 

1/ 

6/ 

30 

PE 

1/ 

6/ 

31 

PE 

1/ 

6/ 

32 

KO 

1/ 

6/ 

33 

PE 

1/ 

6/ 

34 

PE 

1/ 

6/ 

35 

KO 

1/ 

6/ 

36 

PE 

1/ 

6/ 

37 

PE 

1/ 

6/ 

38 

KO 

1/ 

6/ 


1 

1 

1 

6 

1 

1 

0 

1 

1 

2 

1 

0 

1 

22 

0 

12 

100000 

1 

looaoo 

OQE+OO 

22 

200000 


12 


-0 -0 12 


1 

1 

-0 

11 


1 

1 

-0 -0 -0 


. 20 00 00 
-0. -0 
11 


,300000 


.200000 
5. 00E+ 00 
6 

.400000 


. 300000 
1. 00E+01 


.200000 


TOTAL PROS. NOT EQUAL TO 1 


3. G3E+00 
0 
1 
0 
1 
1 

0 
1 
1 

0 

1 
1 

0 
1 
1 

G 

l 


1. 0 0.E+01 

-0 


3.0QE+01 


-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-c 

-0 

-0 


-0 

-0 


-0 -0 -0 


MARINER V ASSEMBLY AND TEST 
SPACECRAFT ASSEMBLY 
INSTALL BAY V SOLAR PANEL 
1.000 1.000 1.000 1.000 

LAMINAR FLOW TENT 

.100 50.000 1.000 60.000 60.000 60.000 60.000 

VISUAL INSP. 

MANIPULATE 3Y HAND 


SOLAR PANEL, BAY V SQX-3755-42Q 

1.000 1.000 


TEN BUMP ID HISTOGRAM 


.150000 

.100000 .050000 

. 040000 

.030000 

. 02000D 

.010000 

0 . 

-0 . -0. 

-0 . 

-0 . 

-0 . 

-0 . 

BASIC 

DISTRIBUTION 





. 150000 

.100000 .050000 

.040000 

.030000 

. 020000 

.010000 

2.00E+01 

5.CQE+01 1.00E+D2 

2.C0E+02 

5.00E+02 

l.QCE+03 

2.00E+03 

PART 

BURDEN SQX3755-420 





.150000 

. G4990C 






1.00E+02 3 

. 00E+02 






PREP OF 

SOLAR PANEL 





-0. 

“C • 

-0. 





10.000 

1.000 






-0. 

-0. 

-0. 

-0. 

-0. 



PART 

1 SURFACE 

ARE AS= 

.100 

0. 

0 . 

0. 

-0 . 

-0 . 

-0. 

-0 . 

-0. 



10.000 

3.000 






-0 . 

-0 . 

-0 . 

-0 . 

-0 . 



PART 

1 SURFACE 

AREAS= 

.100 

0. 

c . 

0. 

-0 . 

-0. 

-0 . 

-0 . 

-0. 



10.000 

2.000 






-0 . 

-0. 

-0. 

-0. 

-0. 



PART 

1 SURFACE 

AREAS= 

.100 

0 • 

c . 

0 . 

-0. 

-C . 

-0. 

-0 . 

-0. 



10.000 

3.000 






-3. 

-0. 

-0. 

-0 . 

-Q . 



PART 

1 SURFACE 

ARt AS= 

.100 

0. 

0 . 

0 . 

-0. 

-0. 

-0. 

-0. 

“G * 



10. COO 

1.000 







Table 4. Typical Output, Tape Alteration Program 



Table 5 Tape 9 Alternation Instructions 


Instruction Format: 

Columns: 1-10* 15 20 21-30* 35 

N K M F J 

*Integers in these fields must be right justified. 

The value of J controls the error message printout; 

J = 0 : No change. (J is initially set to 1.) 

J = 1 : Print fatal errors only. 

J = 2 : Print all errors. 

K determines the type of change to be made: 

K = 0 : Ignore all tape records after record N. 

K = 1 : Replace record N with a new card 

K = 2 : Add M cards after record N. 

K = 3 : Delete M records starting with record N. 

K = 4 : Replace parameter M on record N by F. 

K = 5 : Multiply parameter M on record N by F. 

K = 6 : Begin replacing parameter M by F starting with record N. 

K = 7 : Begin multiplying parameter M by F starting with record N. 

K = 8 : Stop altering parameter M. (Record N is the first 

record not to be altered. ) 

NOTES: 

1) 

2 ) 

3) 

4) 

5) 

6 ) 


Instructions must be in order of record number (N); 

No changes are made in parameters on cards read in; 

An instruction with K=0, 1, or 2 must be the last for 
that record number; 

Cards to be added must be placed directly behind 
the instruction; 

No re-start control card may appear in data on 
Tape 9; restarts must be on cards only. 

To read from cards only (no previous data tape), use an 
instruction card with N-0 r K=2, and M=99999 (or any 
number larger than the number of cards to be read). 

A K=0 instruction card should be used whenever only part 

of Tane 12 is to be used; merelv setting- KK=0 on a 
control card is not sufficient. 


7) 
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As the data from Ta5f©12and from cards is combined to form the new data 
tape, the records are checked for errors such as missing data or data out of 
range. This data checking also maintains the current surface areas of each 
affected part, and these are printed each time the corresponding part is 
mentioned. Whenever an error occurs, an. error message is printed, and a 
count is kept of the number of errors of each type. (See Table 6 for the 
listing of error types. ) Errors are either fatal (e.g., card out of sequence) 
or non-fatal (e.g., area of contact is greater than surface area). Fatal 
errors are those that will cause an error exit from the burden prediction pro- 
gram, whereas non-fatal errors will not stop the prediction program but do 
indicate a problem in the data. When the tape alteration program is combined 
with the prediction program, any fatal errors in the data will terminate the 
job before running the prediction program. 

An advantage of the uniform tape format is that an out-of-sequence record 
will seldom cause termination of the processing; a new Ta_pe 9 will be written 
in spite of errors. Two exceptions are the non-uniform histogram record and 
any data read in from cards (which do not have uniform formats and are not 
identified as to record type). 

Since some data tapes exist that were written in a different format 
(essentially card images), an extra subroutine is provided for converting such 
tapes to the new format. This is necessary since the burden prediction pro- 
gram has been modified to accept the new tape record formats and can no longer 
use the old data tapes. 
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Table 6 Data Check Error Types 

— 1 

1) Index out of range; j 

I 

2) Prerequisite missing; 

3) Necessary data missing; J 

k) Too many cards in set; 

5) Index used previously; 

6) Data out of range; 

7) Card out of sequence; 

8) Affected part area zero; 

j 

9) Contact area greater than surface j 

area j 

10) Remaining area negative; 

j 

11) Total probability not equal to 1,, i 

i 

Note: Errors 1 through 7 are fatal; ; 

| 

I 

Errors 8 through 11 are warnings. 

I 

___ _ L 
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3. CONTROL OF THE NUMBER OF HISTOGRAM INTERVALS 

The maximum number of intervals to be used in histogram calculations is 
adjustable in the range from zero (i.e., a constant) through ten. This is 
accomplished by specifying as 15 on any control card (CC) the identification 
number of a histogram having the desired number of intervals. This "model" 
histogram also provides the probability levels to which burden histograms are 
adjusted. (See Ref. 4 , p 35-56 for discussion of the need for this.) 

Subroutine HAS, the Histogram Adjusting Subroutine, was written to alter 
the number of intervals and the probability levels of any specified histogram 
to match the number of intervals in the model histogram. This is done by con- 
structing the cumulative probability function F(x) of the histogram to be 
altered and determining by interpolation the value of x for which F(x) equals 
each desired probability level. (A numerical example is given in Ref. 4, 
p 50-51*) This routine was originally programmed as a part of HCS, the Histo- 
gram Combining Subroutine, and was separated and modified to serve in the 
present task. 

Although the burden prediction program will run without having any 
model histogram specified, it is best to provide one. Assuming that this has 
been done, all input histograms are checked to see if they have too many inter- 
vals, and if so, they are adjusted by subroutine HAS. One exception to this 
is that the general purpose distributions (as opposed to histograms repre- 
senting burdens on parts) are preserved in their original form in separate 
storage locations before they are adjusted by HAS. This permits these histo- 
grams, which must have identification numbers in the range 1 through 20, to be 
reconstituted without loss of accuracy if a larger number of histogram inter- 
vals is specified later in the run. 


Whenever a new model histogram is specified, the first 20 histograms are 
retrieved from their storage locations in their original forms and adjusted 
by subroutine HAS. If the new model histogram has an equal or larger number 
of intervals than the previous one, this is all that is done. However if the 
new model histogram has fewer intervals, all stored histograms are checked and 
are adjusted if they have too many intervals. 

Note that if 15 is changed on a control card to specify a new model 
histogram and if that histogram has not yet been read in, it must be the first 
histogram to be read in since it must be used to adjust all subsequent 
distributions. 

k. PRINT OUT OPTIONS 

The essence of this task was a reformatting of the printout of the burden 
prediction program to conserve space, A notable modifi-atioc. is that all histo 
grams are now printed on three lines regardless of the number of intervale they 
contain. Many histograms that formerly were printed are now available only on 
TAPE 11 since their information was rarely used. Some compression of printout 
such as the Task Summary was made, but these still remain rather spread out for 
ease of reading. 

Optional printout at the operation and part level is available, and as 

somewhat crowded when printed, being confined to one line per operation or part 
(See Table 7 for a sample of this printout.) The printout option is controlled 
by setting parameter II on any control card (see page 28) and may be changed as 


desired. 



STAGE 7 S/C ASSEMBLY +■ TEST , PASAOENA 
TASK 99 ALIGNMENT CHECKS 


SUBTASK 1, INSERT GUIDE PIN 

OPERATION 1, ENVIRONMENT 1, FROM 264.91 TO 274.91 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, ARE A = .100, TOOL BURDEN= 

OPERATION 1, ENVIRONMENT 1, FROM 274.91 TO 234.91 HOURS, 3 MEN 

PART 5 SL,X-3762-599, SURFACE 1, AR£A= .500, TOCL BURDEN* 

PART 1 Sax-3 755-420 , SURFACE 1, AREA* .100, TOCL BURDEN* 

OPERATION 1, ENVIRONMENT 1, FROM 284.91 TO 294.91 HOURS, 2 MEN 

PART 1 sax-3755-420, SURFACE 1, AREA* .100, TOOL BUROEN= 

OPERATION 1, ENVIRONMENT 1, FROM 294.91 TO 304.91 HOURS, 3 MEN 

PART 1 SQX-3 755-420 , SURFACE 1, ARE A* .100, TOCL BURDEN* 

OFERATION 1, ENVIRONMENT 1, FROM 304.91 TO 314.91 HOURS, 1 MEN 

PART 5 SQX-3762— 599, SURFACE 1, AREA* .500, TOCL BURDEN* 

PART 1 sax-3755-420, SURFACE 1, AR£A= .100, TOOL BURDEN* 

SUBTASK 2, INSTALL GAGE COLLAR 

OPERATION 1, ENVIRONMENT 1, FROM 264.91 TO 274.91 HOURS, 1 MEN 

PART 1 SaX-3755-420, SURFACt 1, AREA* .100, TOOL 8URQEN= 

SUBTASK 3, INST LATCH ASSEMBLY 

OPERATION 1, ENVIRONMENT 1, FROM 404.93 TO 414.93 HOURS, 2 MEN 

PART 1 SUX-3755-42Q , SURFACE 1, AREA* .100, TOOL BURDEN= 

PART 5 Sax-3 762-599 , SURFACE 1, ARE A= . 500 , TOOL BURDEN* 

OPERATION 1, ENVIRONMENT 1, FROM *14.93 TO 424.93 HOURS, 3 MEN 

PART 1 SaX-3755-420, SURFACE 1, AR£A= .100, TOOL BURDEN* 

PART 5 Sax-3 762 -599 , SURFACE 1, AREA= . 500 , TOOL BURDEN= 

OPERATION 1, ENVIRONMENT 1, FROM 424.93 TO 434.93 HOURS, 3 MEN 

PART 1 SaX-3755-420 , SURFACE 1, AREA= .100, TOCL BURDEN= 

OPERATION 1, ENVIRONMENT 1, FROM 434.93 TO 444.93 HOURS, 3 MEN 

PART 1 SaX-3755-420 , SURFACE 1, AREA= .100, TOOL BURDEN* 

SU8TASK 5, ADJUST LATCH ASSY 

OPERATION 1, ENVIRONMENT 1, FROM 498.84 TO 508.84 HOURS, 3 MEN 

PART 1 SaX-3755-420, SURFACE 1, AREA= .100, TOOL BUROEN= 

OPERATION 1, ENVIRONMENT 1, FROM 508.84 TO 518.84 HOURS, 1 MEN 

PART 5 SQX- 3762-599, SURFACE 1, ARE A= .500, TOOL BURDEN= 

PART 1 Sax-3755-420, SURFACE 1, ARE A= .100, TOOL BURDEN= 

SUBTASK 9, POSIT LATCH ASSY/JIG 

OPERATION 1, ENVIRONMENT 1, FROM 444.93 TO 454.93 HOURS, 1 MEN 

PART 1 sax-3758-420, SURFACE 1, AREA= .100, TOOL BUROEN= 

OPERATION 1, ENVIRONMENT 1, FROM 45*. 93 TO 464.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, ARE A= .100, TOOL BURDEN* 

OPERATION 1, ENVIRONMENT 1, FROM 464.93 TO 474.93 HOURS, 1 MEN 

PART 1 sax-3755-420, SURFACE 1, ARE A= .100, TOOL BUROEN= 

OPERATION 1, ENVIRONMENT 1, FROM 474.93 TO 484.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, AREA* .100, TOOL BURDEN= 

SUB T ASK 19, REM OAMPER FROM JIG 

OPERATION 1, ENVIRONMENT 1, FROM 484.93 TO 494.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, AREA* .100, TOOL BU RDEN = 

OPERATION 1, ENVIRONMENT 1, FROM 494.93 TO 504.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, AREA* .100, TOOL BURQEN= 

OPERATION 1, ENVIRONMENT 1, FROM 504.93 TO 514.93 HOURS, 1 MEN 

PART 1 Sax-3755-420 , SURFACE 1, AR£A= .100, TOOL BURD£N= 

OPERATION 1, ENVIRONMENT 1, FROM 514.93 TO 524.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, AREA= .100, TOCL BUROEN= 

OPERATION 1, ENVIRONMENT 1, FROM 524.93 TO 534.93 HOURS, 1 MEN 

PART 1 Sax-3755-420, SURFACE 1, AREA= .100, TOCL BURDEN= 
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Bo Description of Subroutines 

The following subroutines are described in three groups; 

1) Tape alteration program; 

2) Burden prediction program; 

3) Detailed printout program. 


Each group is an independent program and can be run separately if desired, 
using the first subroutine listed as the main program. All are written in 
FORTRAN IV and are as machine -independent as possible. Communication between 
these programs is by means of disk or tape files: 


Tape 2; 
Tape 5s 
Tape 6; 
Tape 9s 

Tape 11: 

Tape 12: 


Binary restart dumps (read and write); 

Card reader; 

Printer; 

Binary data written by the tape alteration program 
for the burden prediction program (read and write), 
Binary data saved from burden prediction for detailed 
printout as desired (read and write); 

Binary data originally written by the tape alteration 
program as Tape 9 hut being modified by that program to 
produce an updated Tape 9= If an older data tape is to 
be converted to the new format, it is also designated 
Tape 12. Tape 12 is read only. 


Whether these files are on disk or tape depends on whether a tape has been re- 
quested or not. If permanent records are desired or if the three programs are 
to be run separately, tapes should be used since disk files are erased at the 


end of each job. 
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1. Tape Alteration Program 

SUBROUTINE TAS , the Tape Alteration Subroutine. This subroutine reads tape 
alteration instruction cards, data cards (if any), and Tape 12 (if used). Tape 
records are read, altered by subroutine RAS, and checked and written on Tape. 9 
by subroutine RG.ifS until the record whose number natches that on the current 
alteration instruction is encountered. At this time the specified change is 
made and the next instruction is read. This continues until a higher numbered 
record is specified in the instruction, at which point the program goes back to 
reading, checking, and writing the tape records. An instruction to start or 
stop changing a parameter value does not accomplish this directly but sets 
appropriate flags in subroutine RAS. See Table 5 for instruction types and 
format. (See Fig. 1 for TAS macro logic.) 

SUBROUTINE CRS , the Card Reading Subroutine. This subroutine reads one card 
under the appropriate FORMAT statement, except for distribution cards (DQ) , of 
which two or four are read, depending on whether the number of intervals is 
less or greater than seven. 

SUBROUTINE TRS , the Tape Reading Subroutine. This subroutine reads one logical 
record from Tape 12 . This record may contain either 14 or 23 words as described 
in Section TT.A.2. 

SUBROUTINE RAS , the Record Alteration Subroutine. This subroutine, by means of 
arrays specifying for each record type the parameters to be changed and the 
change to be made, makes all parameter changes specified by any previous tape 
alteration instruction. These changes continue to be made until another in- 


struction cancels them. 
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SUBROUTINE ROWS, the Record Checking and Writing Subroutine. This subroutine 
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makes all checks that the old Data Check Program used to make (see Table 6 for 
the error types). The next record type to be read (either from cards or .Tape 12) 
is determined, and the record is then written both on the printer and on Tape 9 . 

SUBROUTINE EPS , the Error Printing Subroutine. This subroutine prints the 
appropriate error message and maintains a count of the number of errors of each 
type. At the end of the run, this subroutine prints the error totals and causes 
a program exit if any fatal errors occurred. 

SUBROUTINE PTS , the Page Titling Subroutine. This subroutine prints parameter 
numbers along the top of each page for ease in identifying parameters to be 
changed . 

2. Burden Prediction Program 

SUBROUTINE BPS , the Burden Prediction Subroutine. This subroutine reads control 
cards (CC) and performs the actions required by them. Most burden calculations 
are performed in subroutine MBS, but burden differences and zone reorganization 
is done in BPS. BPS also performs the restart functions. (See Fig. 2) 

SUBROUTINE MBS , the Microbial Buildup Subroutine. This subroutine reads input 
data for task, subtasks, and operations and computes the microbial burden on 
each affected surface as a result of fallout, contact or decontamination. See 
Ref. 2, p 28 for a discussion of this subroutine. (See Fig. 3) 

SUBROUTINE HCS , the Histogram Combining Subroutine. This subroutine calculates 
the probability distribution (in the form of a histogram) of the sum, difference, 
product, quotient, or maximum value of two histograms. For all but the maximum 
value, this subroutine requires both input histograms to have their correct mean 
values in DR(N,l) and will adjust the resulting histogram so that its mean has 
the correct value. HCS always returns the mean value of the output histogram in 




Burden 



iction Program Macro Logic 
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DR(N,l). See Ref. 2, p 15 for a discussion of this subroutine. 


SUBROUTINE HAS , the Histogram Adjusting Subroutine. This subroutine adjusts 
the number of intervals and the probability levels of a given histogram to 
match those of the "model" histogram. (See Section II. A. 3*) 

SUBROUTINE HES , the Histogram Equating Subroutine. This subroutine sets the 
values of one histogram equal to those of another. 

SUBROUTINE HMS , the Histogram Multiplying Subroutine. This subroutine multi- 
plies the abscissae and the mean of a given histogram by a constant. 

SUBROUTINE HWS , the Histogram Writing Subroutine. This subroutine prints the 
mean, probability values, and abscissae of a specified histogram. 

SUBROUTINE MAS , the Microbial Accretion Subroutine. This subroutine updates 
the burden on parts that have not been operated on for some time. 

FUNCTION AVF , the Average Value Function. This function computes the mean of 
a given histogram. 

FUNCTION ZF , This function computes the sum, difference, product, quotient, 
or maximum value of two given numbers (Used by HCS). 

3« Detailed Printout Program 

SUBROUTINE DPS , the Detailed Printout Subroutine. This subroutine reads 
detailed printout instruction cards and controls the printing. Since the 
normal condition is not to print, an instruction is necessary for each task 
or subtask to be printed. (See Table 2 for detailed printout instructions. 'I 



SUBROUTINE PLS , the Parameter labeling Subroutine. This subroutine prints 
the parameters stored on Tape 11 with proper identifying information. 

SUBROUTINE HPS , the Histogram Printing Subroutine. This subroutine prints 
the histogram stored on Tape n. 

4 - Subroutines No Longer Used 

1) The Data Check Program and its subroutine ECS, the Error Count 
Subroutine, are no longer used because their function is included 
in the Tape Alteration Program. 

2) Subroutine HNS, the Histogram Numbering Subroutine, was incorporated 
into subroutine MBS since it was only called twice. 

3) The old restart subroutine, MBRS, is not needed because of the new 
restart capability. 

4) The Mean Value Subroutine, MVS, was replaced by the Average Value 
Function, AVF, to allow greater flexibility in its use. 

5) The Mean and Variance Subroutine, MAVS, is not used because AVF 


now calculates mean values and the variance was not used 
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Cs INPUT DATA 

Table 8 lists all types of input data cards and shows their arrangement. 
Alternate branches are indicated by the flow lines on the left with the index 
values that cause the program to take that branch. The brackets on the right 
indicate cards or groups of cards that are repeated until a blank card occurs. 
The contents and format for each data card are described below. The data are 
identified by the variable names used in the computer program (e.g. , KK, IK). 
Their significance is discussed in Reference 2, page 28; and the notation 
used for these variables in the reference is given in parenthesis in the 
following list. 

Note that each integer constant must be punched next to the right hand 
side of its field (right adjusted). 

CC CONTROL CARD 

KK, IK, II, 12, 13, 14, 15, 16 
FORMAT (815) 

KK=0 causes the program to terminate without reading any 
additional data. 

KK=1 indicates that a RUN is to be initiated and calls for RUN, 

STAGE, and TASK DESCRIPTION CARDS and TASK data inputs. 

KK=2 indicates that a STAGE is to be initiated and calls for 
STAGE and TASK DESCRIPTION CARDS using TASK data inputs. 


ED ENVIRONMENTAL DESCRIPTION CARD 

EQ ENVIRONMENTAL QUANTITIES CARD., 


OD OPERATION DESCRIPTION CARD 

OQ OPERATION QUANTITIES CARD_ 


PD PART DESCRIPTION CARD 

PQ PART QUANTITIES CARD_ 


DD DISTRIBUTION DESCRIPTION CARD 

DQ DISTRIBUTION QUANTITIES CARD (1 or 2 ) 







KK=3 indicates that a TASK is to be initiated and calls for 
a TASK DESCRIPTION CARD and TASK data inputs. 

For KK=1, 2, 3, IK determines the point at which TASK inputs 
start. These inputs and the values of IK for which they are 
read are: 

IK=1 ENVIRONMENTS (EM and all subsequent data) 

IK=2 OPERATIONS (OD and all subsequent data) 

IK=3 PARTS (PD and all subsequent data) 

I K=b DISTRIBUTIONS (DD and all subsequent data) 

IK=5 SUBTASKS (KD and all subsequent data) 

II, 12, 15 , and l6 are used as follows (13 and l4 are not used) 

11 controls the level of detail in printout: 

=0 No change 

= 1 Print subtask information 

= 2 Start operation level printout 

= 3 Start part level printout 

12 controls the restart card dumps: 

=0 No change 

= 1 Restart dump at end of rim only 

= 2 Restart dump at end of each task 

=3 No restart dumps. 

15, if non-zero, is the number of the distribution to be used 
for adjusting the probability levels (See Ref. 4, p 50-51) • 


16 , if non-zero, is the new task environment number. l6 must be 
specified if any fallout is to be determined, but it need be specified 
only on the first control cards 

KK=4 calls for determining and printing the burden difference between 
the burdens stored in special storage distributions II and 12. (These 
distributions have a numbering sequence separate from the distributions 
specified by DD cards.) Storage in these distributions is accomplished 
by specifying the distribution number (1 through 20) on the TASK 
DESCRIPTION CARD (TD). 

KK-5 calls for the printing of areas and microbial burdens by zones. 
IK=0 indicates that the zones are to be defined by 
a deck of ZONE DEFINITION CARDS. 

IK=1 indicates that zones are the same as parts and 

that no ZONE DEFINITION CARDS are necessary. This 
card is to be followed by another CONTROL CARD. 

KK=6 calls for a restart and requires a previously written restart 
dump tape to be loaded on Tape 2. 

A restart contol card should be the first card in the data (card 
or tape) since any previous data is incorporated in the restart dump. 
When encountered by the burden prediction program, a restart control 
card will cause a restart tape to be read. 
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ED RUN DESCRIPTION CARD 

KR, RUN 

FORMAT (15, IX, 7Ab) 

A RUN DESCRIPTION CARD is the first card to follow the 
CONTROL CARD initiating a run. .KR is the run number and RUN 
is a 42 alphameric character description of the run. 

SD STAGE DESCRIPTION CARD 

KS, STG 

FORMAT (15, IX, 7A6) 

A STAGE DESCRIPTION CARD follows a RUN DESCRIPTION CARD or 
a CONTROL CARD initiating a Stage. KS, which must be in the 
range 1 to 20 is the stage number and STG is a 42 alphameric 
character description of the Stage. 

TD TASK DESCRIPTION CARD 

KT, TGK, LI, L2, L3, L4 
FORMAT (15, IX, 7A6, 2X, 415) 

A TASK DESCRIPTION CARD follows a STAGS DESCRIPTION CARD 
or a CONTROL CARD initiating a task. KT, which must be in the 
range 1 to 100, is the task number and TSK a 42 alphameric charac- 
ter description of the task. LI and L2 indicate the finish time 
special storage distributions of prerequisite tasks which are 
to determine the start time of this task. If only one prerequi- 
site (Ll) is specified, the new start time is set equal to the 
time stored in Ll. If no prerequisite is specified, the start 
time is set equal to zero. L3 is the index of the finish time 



special storage distribution if this task is to be prerequisite 
for any later task. L4, if specified, causes the total microbial 
burden at the end of the task to be stored in special storage 
distribution L4 and to be available for later computation of a 
burden difference. IK on the preceding CONTROL CARD indicates 
the card type to follow. LI, L2, L3, and L4 are in the range 1 
to 20. 

EM ENVIRONMENTAL AVERAGE LIFETIME MODIFIERS CARD 
AES(l), AES ( 2 ) , AES (3) , ASS(4) 

FORMAT (4E10.3) 

The ENVIRONMENTAL AVERAGE LIFETIME MODIFIERS appear on the 
first card of the environments inputs. AES(J) is the factor which 
modifies the quantity AET(J) on card EQ. 

ED ENVIRONMENTAL DESCRIPTION CARD 
N, DSC 

FORMAT (15, IX, 4a6) 

Each of the (up to 10) environments requires two input cards. 
The first of these is the ENVIRONMENT DESCRIPTION CARD which gives 
the index N and the 24 alphmeric character description DSC. The 
last such description card is to be blank to indicate that all 
environment inputs have been read. All but the last (blank) 
ENVIRONMENT DESCRIPTION CARD are to be followed by an ENVIRONMENT 
QUANTITIES CARD. 

EQ ENVIRONMENT QUANTITIES CARD 

IEC(N), IET(N) , AEC(N), AET(N), AED(N), AEF(N,l), AEF(N,2), 
AEF(N,3), AEF(N,4) 

FORMAT (215, 7E10.3) 
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IEC(N),AEC(N) describe the background number of organisms/cu ft (C) 
in environment I!. I EC is a "shape" distribution index and AEC is 
the mean of the distribution for this application. 

IET(N) , AET(.N) are the average lifetime distribution shape and mean (V). 

AED(N) is the mean number of organisms per cubic foot of air per man 
in this environment (d). 

ASF(N) is th. fallout velocity mean value (f). 

OD OPERATION DESCRIPTION CARD 
N, DSC 

FORMAT (15, IX, 4a6) 

Each of the (up to 20 ) cataloged operations requires two 
input cards. The first of these is the OPERATION DESCRIPTION CARD 
which gives the index N and the 24 alphameric character description 
DSC. The last of these cards is to be blank to indicate that all 
operations inputs have been read. All but the last (blank) card 
are to be followed by an OPERATION QUANTITIES CARD. 

OQ: OPERATION QUANTITIES CARD 

IOT(N), IOQ(N), IOC(N) 

FORMAT (315) 

IOT(N) is the distribution shape for the time required to perform 
operation N. 

IOQ(N) is the distribution shape for the number of organisms per 
cubic foot per man. 

IOC(N) is the distribution shape for the number or organisms on 
hand or tool (contact contamination). 
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PD: PART DESCRIPTION CARD 

N, DSC 

FORMAT (15, IX, 6A6) 

Each of the (up to 120 ) initial parts requires two input cards„ 
The first is the PART DESCRIPTION CARD giving the part number N and 
the 36 alphameric character description DSC. The last 12 characters 
(in columns 31 through 42) are retained as a permanent serial de- 
scription of the part. The last of these cards is to be blank to 
indicate that there are no more parts inputs. All but the blank 
card are followed by a PART QUANTITIES CARD. 

PQ: PART QUANTITIES CARD 

IAB(N,l) , IAB(N,2), IAB(N,3), IAB(N,4), IAB(N,5), IAB(N,6). 

AAG(N) , AAS(N) 

FORMAT (615, 2E10.3) 

IAB(l) is the index of the distribution of microbial burden on 
the top surface of the part. The area of this and the 
other surfaces is given in the coefficient associated with 
the distribution (see the explanation for DQ below). 

IAB(2) is the index of the distribution for area and burden of 
all other exterior surfaces. 

IAB(3) is the index of the distribution for area and burden of 
the mated surface. 

IAB(4) is the index of the distribution for area and burden of 


the occluded surface. 
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I\B(5) is the environment of the part used in updating microbial 

accretion for parts that have not been operated on for some 
time . 

IAB(6) is the distribution shape for part contact retention (Si). 

A--iG is the surface retention factor for fallout contamination (g). 

AAS is the mean part contact retention (ol). 

DD: DISTRIBUTION DESCRIPTION CARD 

N, DSC, L, M 

FORMAT (15, IX, 4a6, 215) 

Arrays are defined for storing up 500 histograms, each with a maximum 
of 11 absissa values and their associated probabilities. Since the program 
automatically indexes new histograms that it must store, only a fraction of 
these 500 should be specified at input time. The DISTRIBUTION DESCRIPTION 
CARD gives the distribution number N and 24 alphameric character description 
DSC. The last of these cards is to be blank to indicate that there are 
no more distributions inputs. 

The index L indicates whether the range values and associated proba- 
bilities are to be read directly from the DISTRIBUTION QUANTITIES CARDS or 
calculated from the range minimum, mode and maximum values. 

M is the number of values (to be in the range 1 to 11) If M=l, the 
distribution is a constant with value XR(N,l). 

DQ: DISTRIBUTION QUANTITIES CARDS 

For L = 1: DR(N,J), J=1,M 

XR(N,J),J=1,M 

For L = 2: DR(N ,1) ,X1 ,X2 ,X3 

FORMAT (8E10o3) 


For L = 1, there are two sets of data DR(N,J) and XR(N,J). DR(N,l) 
is equal to the area for distributions describing the microbial burden on 
particular surfaces. For other distributions it is left blank and the 
program automatically supplies the mean of the distribution. DR(2) to 
DR(M) are the probabilities associated with the intervals XR(l) to XR( 2), 
XR(2) to. XR(3),...XR(M-l) to XR(M). If M is larger than 8, two cards are 
required for each set of array values. 

For L = 2, DR(l) has the same significance as for L=l. The range 
intervals of the distribution and the associate probabilies are deter- 
mined by the program to describe, with a 5 value histogram, a triangular 
distribution with the range XI to X3 and mode (probability peak) at X2„ 

The program automatically generates and indexes (index = KO), a zero 
value histogram. 

SUBTASK DESCRIPTION CARD 

N, DSC, Nl, N2 

FORMAT (15, IX, 4a6, 215) 

The input data for each of the (up to 20) subtasks per task is 
headed with a SUBTASK DESCRIPTION CARD that gives the subtasks number 
N, the 24 alphameric character description DSC, and the indices of two 
prerequisite subtasks, Nl and N2. If Nl and N2 are blank, the subtask 
start time is taken to be the same as the task start time. If there is 
only one prerequisite subtasks, it is to be indicated by Nl since, if Nl 
is blank, N2 is not checked. Both subtasks must have been previously 
considered by the program and had their finish times computed. 
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KG: S'TBTASK CHANGjS CARD 

K, IE(1), IR(2), Ifi(3), IR(4) , IR(5), AR(l) , AR(2), AR ( 3 ) 

FORMAT (615, 3E10.3) 

The SUBTASK CHANGE CARD is used to define any changes in area, 
environment, or retention factor values for a part. The amount of area 
AR(1), is subtracted from part IR(3)» surface IR(4), and added to part 
1 •: ( 1 ) , surface IR(2). 

If IR(5) , is a positive integer, it is the new environment number 
assigned to part IR(l). 

If AR(2) is a non-zero number, it is the new surface retention 
factor (g) for fallout for part IR(l). 

If AR(3) is a non-zero number, it is the new mean part contact 
retention factor (Si) for contact for part IR(l). 

K is the change index and is to be an arbitrary positive integer 
except where it is to indicate there are no more changes, in which case 
it is left blank. 

KO: SUBTASK OPERATION CARD 

10, IKE AKT, AKQ 
FORMAT (215, 2KL0.3) 

Each subtask consists of an arbitrary number of operations. 

These operations are partially described in the catalog of operations 
(see inputs 0D and OQ described above) but have additional parameters 
to describe the resulting microbial burden changes for each particu- 
lar performance of the operation. 

10 is the operation number. If 10 is zero or blank, there are no more 
operations and the program proceeds to a KD card. If 10 is negative, 
a decontamination is indicated for which the input quantities have a 
a different significance. 



For decontamination (ICXO): 


IKE is not significant and can be left blank. 

AKT is the operation time. 

For fallout and contact (10 >0): 

IKE is the operation environment unless it is left blank, 
in which case the operation environment is the same as 
for the previous operation. If IKE is blank for all 
operations in a subtask, the operation environments are 
all the same as the environment number l6 given in card CC. 

AKT is the mean operation time and is used to modify the distri- 
bution I0T(I0) (see card OQ) to determine the operation time 
distribution. 

AKQ is the number of men close to the hardware during this 
operation. 

PART EFFECT CARD 

IP, LS, APD, APC , APS, APA(l), APA(2) or IP, LK, AR(1), AB(2), AR(3), 

AR(4) 

FORMAT (215, 5E10o3) 

IP is the number of the part affected. If IP is blank or zero, 
it indicates that all affected parts have been considered and 
the program is to proceed to a consideration of the next 
operation. 

For fallout and contact operations: 

LS is the tool retention distribution for contact. 

APD is not used. 


APC is the mean number of organisms on hand or tool (contact 

contamination) . 
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APS is the retention factor for the tool in determining 
contamination by contact (S2). 

APA(l) is the contact area for the top surface, APA(2) the contact 
area for other exterior surfaces. It is assumed that occluded 
and mated surfaces are not contacted. 

For decontamination operations: 

LK is the distribution shape of the fraction of organisms removed. 

AR(J) is the mean fraction of organisms removed from surface J. 

ZD: ZONE DESCRIPTION CARD 

IZ, DSC 

FORMAT (15, IX, 4A6) 

Defintion of each thermal zone in terms of the fractions of 
surfaces that comprise it requires a ZONE DEFINITION CARD for the 
zone and one ZONE COMPOSITION CARD for each part involved. The zone is 
described in terms of a number, IZ, and a 2k alphameric character de- 
scription DSC. If IZ i6 blank or zero it indicates that all zones 
have been defined and that the program is to proceed to the next 
CONTROL CARD. 

ZC: ZONE COMPOSITION CARD 

IP, FP(l) , FP(2), FP(3) » FP(4) 

FORMAT (15, 5X, 4F10.7) 

A separation ZONE COMPOSITION CARD is used for each part IP 
that is wholly or partially included in zone IZ. The quantities FP(J) 
indicate the decimal fractions of surface J (top, exterior, mated, and 
interior) to be included in zone IZ. The last ZONE COMPOSITION CARD 
for each zone must be blank; this is followed by the next ZD card. 
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D. JOB CONTROL AN D DECK ARRANGEMENT 

Table 9 lists a typical input deck for the Microbial Burden Prediction 
Program, the control cards are for the Univac 1108 computer. 

It is assumed that a previous run has written a data tape (6524) and that 
it is desired to use this tape, add some more data from cards, and compute 
a burden prediction. Detailed information from the burden calculations is 
to be written on Tape 6523? and is to be printed for the last task computed 
(Stage 5i Task 7)* A new data tape (6522) is also to be written. 

Tape 6522 is designated Tape 9 and will replace 6524, which is desig- 
nated Tape 12. (Tape 6524 was designated Tape 9 when it was written during 
the previous run.) Tape 6523 is designated Tape 11, and will contain de- 
tailed information that can be selected for printing at some later date. 

The main program contains only four call statements; it calls each of 
the major subprograms (Tape Alteration, Burden Prediction, and Detailed 
Printout) in turn and exits. The remainder of the source deck is as listed 
in the Appendix. 

The first data card is read by the Tape Alteration Subroutine (TAS). 

It indicates that parameter number 3 on record number 103 is to be changed 
to a 7- The result of course depends on what type of card record 103 
represents; it could, for example, be a change in the number of histogram 
intervals to be used in the computations. The second data card is also read 
by TAS; it causes (up to) 999 cards to be read in after record number 1755* 
These records and the records obtained from Tape 12 are written on Tape 9 to 
form an updated data tape. Note that it is not necessary to count the cards 
to be added as long as they are added to the end of the deck. Any number 



Table 9« arrangement of data deck 


&RUN?/TPC LBB? 06024S? 07695 » 05» 199/000 , MBP L BUSCH 125/45 

aMSG MOUNT TAPE 6 522 TO WRITE LBB 
AMSG MOUNT TAPE 6523 TO 'WRITE LBB 
AA56 » T 9?T?6522W ) 

MSG?! UrT»6523W 
AASG^T 12?T?6524R 
afor 

CALL TAS 
CALL BPS 
CALL DPS 

CALL EXIT • 

END 

AFOR 

SUBROUTINE TAS 

— REMAINDER OF SOURCE DECK — — 


AXQT 


4 3 

2 999 


7 


c A ^ 

7 s/C ASSEMBLY 4- TEST p PASADENA 
5 INSTALL BAY VII SOLAR PANEL 
1® A# 1000 e 5 , 

5 HIGH BAY 


100, 


5 SOLAR PANEL? BAY Vi I 
jS 

0 

6 ONE BUMP ID HISTOGRAM 

i • 


5QX“3762~5 t 

1 1, 


5 HIGH BAY SURFACE DlSTR 
1 , 

0 99 1 s 

25 SOLAR PANEL? EXTERNAL 

1 0 1 e 

*99 1. 

0 

5 XNSP OF SOLAR PANEL 


i 

1 


REMAINDER OF DATA DECK 


(TAS instructions 
■ctas instructions 

(CONTROL CARDS 

li 2 


100, 1OQ 0 


(CONTROL CARD! 
(DPS INSTRUCTIONS 
COPS INSTRUCTION) 

aPMD » AE 
aFIN 


15 7 

0 


NO" 


rHA 


.presents a master 


(MULTIPLE 7* 


HJN 
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larger than the number of cards to be read can be used for the parameter M on 
the TAS instruction; reading of these cards will be terminated after the last 
control card since this control card indicates the end of the data deck. (By 
inference, record number 1756 on Tope 12 was a control card indicating the end 
of the previous data deck. This record does not appear on Tape 9> being re- 
placed in the record sequence with a different control card.) 

Immediately after the last control card (the last card read by TAS) are two 
instructions to the Detailed Printout Subroutine (DPS). The first of these in- 
dicates that all detailed information resulting from burden calculations on 
Task 7 of Stage 5 are to be printed. The last DPS instruction amounts to a 
program stop; its omission would cause an abnormal exit since DPS reads another 
instruction as soon as it has completed executing one. 
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III. RE SULTS AND CONCLUSIONS 

In general, all objectives of Phase VII were met; although the results 
in some cases differed from what was expected. 

The restart capability performed satisfactorily; but since no restart in- 
formation is provided for the tape alteration program, a number of false error 
messages (e.g., "prerequisite missing") are generated during data checks,. These 
error messages can of course be ignored but rnay possibly obscure a valid error 
message. A restart capability could be added to the tape alteration program, 
but it would require an extra tape and was therefore decided against. 

The detailed printout capability performs well and is easy to use; it 
promises to be a very useful addition to the program. The ability to dispense 
with tapes for short runs and to read all inputs from cards is a great convenience. 

The tape alteration program does everything that was expected of it. This 
capability should be particularly useful for sensitivity studies. 

Although no extensive comparisons were made, it appears that the running 
time of the burden prediction model does vary as the square of the number of 
histogram intervals used. Calculation of the Mariner Venus 67-2 assembly and 
test sequence (using data prepared in Phase VI) required 20 minutes using 5~ 
interval histograms and seven minutes using three intervals. 

The re-formatted outputs were somewhat disappointing since the print-out 
volume was reduced only by approximately 20%, This is because the addition of 
a complete task summary (burden histograms for each surface of each part) al- 
most offset the reduction obtained by re-formatting. It was concluded that 


such detailed task summaries should be made optional„ 
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No trouble was experienced in getting the computer program to run on the 
Uni vac 1108 at JPL. (It was originally checked out on the CDC 6400 at Martin 
Marietta's Denver facility.) However, some modification of the program was 
necessary to make effective use of the JPL computer. In particular, the restart 
dumps were changed from punched cards to tape. 

Although the purpose of Phase VII was not to produce a burden prediction, 
a prediction for the Mariner Venus 67-2 spacecraft was generated by the demon- 
stration case (Phase VI data). Comparison of the Phase VII prediction mean 
value with that obtained in Phase VI was made at two points during the runs: 



Phase VI Phase VII 

Stage 1, Task 8l 

(End of Pasadena Operations) 

121600. 202000. 

Stage 2, Task 15 
(Burden at launch) 

79700. 207000. 


Most of this difference is attributed to an error that was discovered in the 
burden prediction model during Phase VII (a reference histogram was being 
changed inadvertently). Part of the difference is due to the higher contact 
burden (corresponding to vegetative organisms) use in Phase VII; it was not 
judged worth the expense to make a special computer run in Phase VII merely 

to alter this parameter. The Phase VII burden prediction was compared with 
two estimates obtained by JPL as described in Ref. 6: 

Long-term cupons (Ref. 6, Table k, May 28): 3.2 x 10^ spores 
Swab-rinse (Ref. 6, Table 6, May 28): 3-9 x 10^ spores 

The Phase VII calculation showed that the burden will exceed these values 
with probability 0.17 and 0.10, respectively. 



44 


TV, RECOMMEN DATIONS FOR PROGRAM CONTINUATION 

The work performed during Phase VII was primarily intended to improve 
the flexibility of the Microbial Burden Prediction Model so that sensitivity 
studies and parameter adjustments could be more easily made. Although the 
Mariner Venus 67-2 data prepared in Phase VI could be used for this purpose, 
the Mariner Mars 69 data are now available and incorporate a number of im- 
proved metho's for bioassay and assembly sequence recording. It is therefore 
recommended that a set of inputs be prepared to simulate the Mariner Mars 69 
assembly and test sequence so that sensitivity and parameter adjustment studies 
can be performed on the best available data. 

Since the routine manual effort of card punching, sorting, and collating 
formed a considerable part of the Phase VI work, it is also recommended that 
a computer program be written to perform the routine work associated with in- 
put preparation. Such a program should also include a complete list of parts 
at the level of detail recorded in the assembly log, and should be able to 
maintain a running account of the status (i.e. , installed or not) of each part. 
This computer program should be prepared with the capability of adaptation to 
a real time burden prediction capability for use on future spacecraft assemblies. 

An additional recommendation is to make the bulk of the Task Summary print- 
out optional; that is, omit printing the zone-by-zone burdens unless requested. 

It may be preferable to permit this option at the subtask level as well. 
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APPENDIX 


PROGRAM SOURCE DECK LISTINGS 



A2 


subroutine tas 

C TAPE ALTERATION SUBROUTINE 

C THIS SUBROUTINE READS FROM TAPE12 AND WRITES ON TAPE9 

COMMON/T/LN * LR f IX (6) »XX(22) »M»Nl»N2» I (6»20) »X(22»20) »K(22»20 ) t 
• KIL » NE » KE ( 1 1 ) f I A6 (120 > 6 ) > DR ( 500 ) » XR ( 500 ) , NU » 10 » AR ( 4 ) , KKT ( 20 5 . NP* 

. LH ( 20 ) »KXT A 

600 FORMAT (/13H CHANGE » 1 10 * 215 , F15 .6/ ) 

KXT-0 A 

DO 1000 L=1 f 11 
1000 KE(L)=0 
Ni=0 
N2"0 
K IL-1 
ND=0 
NE=0 
LN=1 

CALL PTS 
DO 1 J=l»22 
DO 1 L“i * 20 

1 K< J»L)=0 

READ ( 5 * 500 ) NN »KK »MM * F » KLL 

IF( (NN.NE.O) .OR. (KK.NE.2) ) GO TO 9 

NP=NP+3 

1F(NP.6E.S5) CALL PTS 
WRITE (6 » 600 ) NN * KK * MM * F 
GO TO 21 

2 READ < 5 » 500 ) NN*KK*MM»F »KLL 
500 FORMAT( 110 » 215 #F10. 0 » 15) 

9 IF(KLL.GT.O) KIL=KLL 

3 CALL TRS 

4 IF(NN-Nl) 6666 * 7 # 5 

5 CALL RAS 
CALL RCwS 

IF(LN.EU.O) RETURN 
GO TO 3 
6666 NP=NP+5 

IF ( NP . GE • 55 ) CALL PTS 
WR I TE ( 6 * 60 1 ) 

601 FORMAT (/34H ***** INSTRUCTION OUT OF SEQUENCE) 

WRITE ( 6 * 600 ) NN»KK»MM»F 

6 READ( 5 » 500 ) NN*KK*WM»F»KLL 
IF(KEL.GT.O) KIL=KLL 

60 TO 4 

7 NP-NP+3 

IF ( NP » GE . 55 ) CALL PTS 
WRI TE (6 * 600 ) NN»KK»MM»F 
IF(KK.GT.O) GO TO 6 
CALL RAS 
CALL RCWS 
CALL EPS ( 0 ) 

RETURN 

8 GO TOdUlsl 02 *103*1 04 *105*1 06 *107 *108)* KK 

C REPLACE RECORD NN WITH A NEW CARD 

101 CALL CRS 
CALL RCWS 

IF(LN.EQ.O) RETURN 
GO TO 2 

C ADD MM CARDS AFTER RECORD NN 

102 CALL RAS 
CALL RCWS 
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IF ( LN * EU » 0 ) RETURN 
21 DO 20 L=1»MM 
CALL CRS 

CALL ROWS ...... 

IF(LN.EQ.O) RETURN 

20 CONTINUE . 

60 TO 2 

C DELETE MM RECORDS STARTING WITH RECORD NN 

103 DO 30 L=1»MM 
CALL'TRS 

30 CONTINUE 
60 TO 6 

C REPLACE PARAMETER MM ON RECORD NN BY F 

104 IF ( LR • EQ • 17 ) 60 To HO 
IF(MM.LE.6) 60 TO 41 
XX(MM-e>)=F 

GO TO 6 
41 IX(MM)=F 
60 TO 6 
40 XX(MM)=F 
60 TO 6 

C MULTIPLY PARAMETER MM ON RECORD NN BY F 

105 IF ( LR « EQ » 17 ) 60 To 50 
IF ( MM « LE « 6 ) 60 TO 6 
XX(MM-6)=XX(MM-6)*F 
60 TO 6 

50 XX-(MM)=XX(MM)*F 
GO TO 6 

C REPLACE PARAMETER MM BY F ON NN AND SUBSEQUENT RECORDS 

106 K ( MM'r LR ) = 1 

60 IF (LR *EQ . 17 ) GO To 66 
IF ( MM • LE * 6 ) 60 TO 64 
X(MM-b»LR)=F 
60 TO 6 
64 I ( MM » LR ) -F 
60 TO 6 
66 X ( MM » LR ) -F 
60 TU 6 

C MULTIPLY PARAMETER MM BY F ON NN AND SUBSEQUENT RECORDS 

10? K(MM?LR)~2 
60 TO bQ 

C STOP ALTERING PARAMETER MM, (RECORD NN IS NOT ALTERED® ) 

106 K ( MM * LR ) -0 
60 TO 6 
END 1 



subroutine trs 

TAPE READING SUBROUTINE (READS BINARY RECORDS FROM TAPE12) 
C0MM0M/T/LN»LR»IX(6) ,XX(22) »M»Nl»N2» I (6 » 20 > * X ( 22 » 20 ) »K < 22 1 20 ) r 
®KIL*NE#KE (11) * IA6< 120 » 6) * DR (500) »XR(500) »ND» I0»AR(4) »KKT<20) *NP* 
•LH(20) »KXT 
N1=N1+1 

IF (LNeEQe 17 ) GO TO 5 
READ(12)LR» IX » (XX(J) »J=1»7) 

RETURN 

5 READ ( 12 ) LR » XX 
RETURN 
END 


At, 

subroutine crs 

C CARD-READING SUBROUTINE 

COMMON/ T/LN » LR » IX ( G ) , XX 1 22) »M»Nl »N2» I (6» 20> » X C22»20> »K(22f20) r 
• K It * NE » KE ( 1 1 ) p IAB(120f6) » OR ( SOO ) » XR t 500 ) f NDf 10 » AR ( 4 ) »KKT<20 ) pNPp 
® EH ( HO ) j KXT A 

LR=LN ' 

GO TO(lr2f3#4»5»6»7i8#9»10»llri2»13»14»15»Ibfi7»18*19>5)»LN 

1 READ (5*501) IX 

501 FQRMAT(4X5»1GX»215) 

RETURN 

2 READ(5p 502> IX (i ) »'( XX ( J > * J=1 » 4 ) » I X ( 2 ) » IX ( 3 ) 

502 FORMAT ( I5r 1X*4A6?2I5) 

RETURN 

3 READ ( 5 * 503 ) IX » ( XX ( J ) » J=1 # 3 ) 

503 FORMAT(615»3E10.3> 

RETURN 

4 READ ( 5 » 505 ) IX (1 ) * IX ( 2 ) » XX ( 1 ) * XX < 2 ) 

RETURN 

5 READ (5*505) IX < 1 ) r IX (2 > r ( XX ( J) * J=1 » 5) 

505 FORMAT (2Iae7ElQ®3) 

RETURN 

6 READ ( 5 f 508 ) IX ( 1 ) » ( XX ( J ) * J=1 » 7) 

RETURN 

7 READ 1 5 » 508 ) IX ( 1 ) »’( XX ( J) » J=1 » 7 ) 

RETURN 

Q READ ( 5 p 508 ) XX ( 1 ) » ( XX ( J ) » J=1 t 7 > » ( I X ( J ) p J=2 » 5 ) * 

500 FORMAT(l5*lX»7Ab»2X»4I5) 

RETURN 

9 READ ( 5 p 509 ) (XX ( J) ♦ J=l»4) 

509 FORMAT (8E10. 3) 

RETURN 

10 READ(5p502) IX(I) » (XX(J) »J=1»4) 

RETURN 

11 RE.AD(5p 505) IX(1) » IX(2) r (XX(J) * J=li7> 

RETURN 

12 READ (5*502) IX (1 ) * ( XX ( J ) * J= 1 r 4 ) 

RETURN 

13 READ(5»501> < IX ( J> # J=1 » 3 ) 

RETURN 

14 READ (5*508) IX ( 1 ) » ( XX ( J ) * J=1 * 6 ) 

RETURN 

15 READ ( 5 o 503) IX * XX ( I ) , XX (2 ) 

RETURN 

16 READ (5*502) IX ( i ) » ( XX < J ) » J=1 1 4 ) » XX ( 3 ) ♦ IX ( 2 ) 

RETURN 

17 READ ( 5 9 509 ) ( XX ( o ) » J=1 » M ) 

MPsM+11 

READ (5 5 509 ) ( XX ( J ) » J=12 * MR > 

RETURN 

18 READ ( 5 p 502 ) IX (1) * (XX (J) » J=1 *4) 

RETURN 

19 R£AQ{5e519) IX < 1 ) # (XX( J) * J=l>4) 

519 F0RMAT(I5»5X#4E10.7) 

RETURN 

END 


SUBROUTINE ras 

RECORD ALTERATION SUBROUTINE 

COMMON/ T/LN #LR » IX (b ) » XX ( 22 ) » M » Nl # N2 # I ( 6 » 20 > # X ( 22 » 20 > t K < 22 » 20 ) » 

• KIL»NE»KE(I1) » I AB ( 120 » 6 > » DR < 500 > # XR ( 500 ) # ND t 10 * AR ( 4 ) #KKT(20) #NP 
•LH(20) #KXT 
IF(LR.EQ.17> GO TO 29 
DO 10 LSI ,6 

IF(K(L»LR)®EQ®1) IX (L)=I (L»LR ) 

10 CONTINUE 

DO 20 L=7#13 
IF (K(L»LR)“I) 20 * 21 #22 

21 XX(L-6)=X(L-b#LR) 

GO TO 20 

22 XX{L-6)=XX(L-6)*X(L-6#LR) 

20 CONTINUE 

RETURN 

29 DO 30 Lsl f 22 

IF (K (L»LR)~1) 30 # 31 » 32 

31 XX(L)=X(L»LR) 

60 TO 30 

32 XX(L)=XX(L)*X(L#LR) 

30 CONTINUE 
RETURN 
END 



Ay 

n 

subroutine rows 

C RECORD CHECKING AND WRITING SUBROUTINE (WRITES BINARY ON TAPE9i 
C LR IS CURRENT RECORD TYPE* LN IS NEXT RECORD TYPE ® 

C0MM0N/T/LN»LR»IXC6) *XX(22) *> Hip Hi * N2 1 1 ( 6 , 20 ) # X ( 22 » 2Q 5 * K ( 22 » 20 1 » 
•KIL»NE»KE<11> »IA8(120»6) * DR (500) eXRCSOO) »ND» I0*ARC4) *KKT<20) *NPe> 

*LH (20 ) » KXT A 

DIMENSION CT (203 

DATA (CT ( I ) # 1=1*20) /2HCC » 2HKD » 2HKC » 2HK0 * 2HPE * 2HRD * 2HSD * 2HTQ » 2HEM * 

. 2HED » 2HEQ » 2H0D » 2HOO * 2HPD * 2HPQ * 2 HDD * 2HDQ * 2H^D * 2HZC * 2HPE/ 

6uO FORMAK I7f 3X#A2» X4»lH/»I3*iH/» 12) 

N2=N2+1 

LX=LN 

NP=NP*1 

IF (LRoEQ. IT ) NP=NP+2 
IF ( NP »GE * 55 ) CALL PTS 

GO TO (10 *20 '30 '40 '50*60 '70 '80 *90 '100*110 '120*130 '140*150 '160 *170* 
a 160' 190*200) »LR 

C CC CARD 

10 WRITE (6 » 600 ) N2#CT(LR> 

WR ITE ( 6 » 601 ) IX 

601 F0RMAT(1H4*23X*6I4) 

KK=1X(1) 

IF (KK • GT • 0 ) GO TO 11 
WR ITE ( 9) LR»XX» ( XX < J) *J=i*7) 

CALL EPS ( 0 ) 

LN~0 

return 

11 IF (KKoLE.5) 60 TO 12 
K X T - 1 ■ 

LN-1 

IFCKK.EQ.6) GO TO 1000 
KXT-0 

CALL EPS ( 1 > 
ln=q 

60 TO 1000 

12 1F(KK~4) 13*14*15 

13 IK"IX(2) 

LN=KK+5 

IFC t IK.GT.Q) .AND. (iKoLE.S) ) So To 1000 
CALL EPS CD 

j K“ r i 

60 TO 1000 

14 XF< ( IX<3) .LE.O) .OR. (IX (3) .GT.20) ) CALL EPS(l) * 

IF ((IX(4> .LE.O) .OR. (IX(4> .GT.20 > > CALL EPSU) 

LN-1 

GO TO 1000 

15 LN=18 

IF (IX (2) .EQ.l) LN-1 
60 TO 1000 

C KD CARD 

20 N=IXC1) 

WRITE (6*600 ) N2'CT(LR) »KS'KT*N 
WRITE ( 6 ? 602 ) ( IX ( J) ' J=1 * 3 ) * (XX ( J) p J"1 » 4 ) 

602 F0RMAT(1H+'23X*3I4'17X'4A6) 

IF ( N * GT • 0 ) GO TO 21 

L N — 1 

SO TO 1000 

21 LN=3 

IF ( N «6T « 20 ) 60 TO 28 
IF C KXT i N ) o EQ , 1 ) CALL EPS C 5 ) 

IF (IX (2) .LE.O) 60 TO 29 


3 ^ 5 » 


M1=XX(2) 

IF (Ml. GT. 20) GO TO' 26 
IF (KKT (MX ) «{£Q .0 ) CALL £PS(2) 

IF ( IX (’3 ) .LE.O) GO TO 29 
M2=IX<3) 

IF(M2*GT.20) GO TO 28 
IF (KKT (M2) .EQ.O) CALL EPS (2 ) 

GO TO 29 
20 CALL EPS ( 1 ) 

29 KKT ( N ) = 1 
GO TO iOOO 

KC CARO 

30 WRITE (6*600 ) N2»CT<LR) »KS»XT»N ' 

WRITE (6 ? 603 ) IX » < XX ( w ) , J=l , 3 ) 

603 FORMAT (1H+*23X» 6IA »3FI0.3) 

IF ( IX ( 1 ) .GleU) GO TO 32 
LN=4 

GO TO 1000 

32 IFCXX(l) .GT.O. ) GO TO 33 
IF(IX(2>> 38*38* 39 

33 X1=IX(2) 

X2=IX(3) 

IF( (Il.LE.O) .OR. (I1.GT.120) ) GO TO 38 
IF ( ( I2.LE.0) .OR. ( I2.6T.4) ) GO TO 38 
IB=IAB(I1»X2) 

I3=IX<4) 

X4=IX<5> o 

IF( U3.LE.0) .OR. < I3.GT.120) ) GO TO 38 
IF( ( I4.LE.0) „0R. ( I4.GT.4) ) GO TO 38 
I A- X AB (13.14) 

XF(IA.GT.O) GO TO 34 
CALL EPS ( 3 ) 

GO TO 39 

34 IF( IB.GT.U) GO TO 35 
NOs NO 41 

IB=ND 

IAB(I1» I2)=IB 

DR C IB ) = 0 . 

35 DR(IA)=0R(IA)-XX<1) 

DR(IB)=OR(IU)+XX(l) 

WR I T£ ( 6 * 633 ) DR(IU) »DR(IA) 

633 FORMAT (1H+»83X»9HNEW ARE A= * F8 . 3 * 5X * 1 OHAREA LEfT=*F8.3) 
IF(DRUA)) 36*37* 39 

36 CALL EPS (10) 

37 I AB ( 13 ? 14 ) so 
GO TO 39 

38 CALL EPS ( 1 ) 

39 LNs3 

GO TO 1000 

KO CARD 

40 WRITE ( 6 * 600 ) N2 » CT C LR ) t KS * KT * N 

WRITE (6» 604) ( IX ( J ) * J=1 * 2 ) * ( XX ( J) * J=1 *2 ) 

604 FORMAT ( 1H+»23X»2I4» 16X*2F10.3) 

10=1 X C 1 ) 

IF C 10 > 49*41*42 
49 LN=20 

GO TO 1000 

4 1 LM~2 

GO TO 1000 

42 IF ( 10.GT «2Q ) CALL EPS(l) 
lF((lX(2).LE.-l).0R.(IX(2).GT.ia)) CALL EPS(l) 



WRITE (6*600) N2»CTILR} t Kb » K T » N 


WRITE (6 t 605) (IX ( J ) # J=1 » 2 ) , ( XX ( J) * J=1 » 5 ) 
FORMAT ( 1H+ c23X»2l4»i6X»5F10.3) 


trn-t Y i 1 T 
i r — I a % X ) 


XF(IP«GTsQ) GO TO 51 
LN=4 

GO TO 1000 

XF( XP.GT.120) GO TO 58 

IF ( IX (2 > oGl <s 500 ) CALL EPS(l) 

AA«0 ® 

DO 54 J=l»4 
AR( J)-0o 
IB-IAB ( IP »<J ) 

IF( IB.LE.O) GO TO 54 
AR ( J ) ~OR ( IB ) 

A A- AA+ AR ( J ) 

CONTINUE 
NP=NP + 1 

WRITE ( 8 <>635 ) IP e AR 

FORMAT I55X»4HPART» X4»15H SURFACE AREAS' 
IF (AAoLEsOs ) CALL EPS(Q) 

00 55 Jsieg 

IF(XX(J+3) .OT.AR(J) > CALL EPS (9) 


EPS (9) 


CONTINUE 


GO TO 59 

CALL EPS ( 1 ) 
LN=5 

eo to iooo 


RD CARD 


: 1 . 7 ) 


WRITE16»600) N2»CT (LR) 

WRITE (6 » 606) IX ( 1 ) # (XX ( J) » J=i . 7) 
FORMAT ( 1H+ 1 23X » 14 » 2bX r 7A6 ) 

DO 61 J- I s 120 
DO 61 L=lf6 
IAB( J? K)=0 
DO 62 J~lf500 
DR( J)=0. 

DO 63 jsl e 20 
LH ( J ) -0 
L.N-7 

60 TO 1000 

SD CARD 

KS=IX(1) 

WRITE (6*600) N2 * CT ( LR ) » KS 
WRITE (6 >607) KS» (XX(j) » J=l» 7) 
FORMAT ( 1H+ * 23X t 14 , 25X * 7A6 ) 

L.N-8 

IF( (KS.QT.O) .AND. (KS.LE.20) ) GO TO 


If { CKSsGToQ) 
CALL EPS Q) 
KS~~1 

60 TO 1000 


1000 


CARD 


:IX(1) 


WRITE (6? 600 ) 


T ( LR ) 


WRITE (6? 608 ) CIXCJ 
FORMAT ( 1H+ 1 23X s> 5X4 
00 SI J=i»20 
J 


?*J=i 


1F((KT*GT.0) » AND« (KT eLEelOO) ) GO TO 82 

CALL EPS ( 1 3 
KT~“1 

A3 t 1-TV/Ol 

IF C Ll « LE « 0 ) GC TO 84 
IFCH.6T.20) GO TO 83 
IF (LH ( LI ) *LE • 0 ) CALL £P3(2> 

L2=IX(3> 

IF (L2 e LE » 0 ) GO TO 84 

IF C L2 ® GT o 20 3 GO TO 83 

IF < LH ( L2 ) *LE • 0 ) CALL EPS (2) 

GO TO 84 

83 CALL EPS < 1 ) 

84 L3-IX ( 4 ) 

IF(L3.LE.0> GO TO 8o 
IF (L3»GT.20 ) GO TO 85 
LH (L3 ) “1 

GO TO 36 

85 CALL EPS ( 1 3 

86 IF ( IX (5) »GT «2U 3 CALL EPS C 1 ) 

GO 10(801 r 802 » 803 e 804 #805) t IK 

801 LN-9 

GO TO 1000 

802 LN=12 

GO TO 1000 
603 LN~14 

GO TO 1000 

804 LN=i6 

GO TO 1000 

805 LN=2 

GO TO 1000 

C EM CARD 

90 WR ITE ( 6 f 600 3 N2 »CT (lR) »KS?KT 
WRITE (6 ? 609 ) ( XX ( o ) , j=l » 4 ) 

609 FORMAT (1H+»47X»4F10.3) 

LN=10 

GO TO 1000 

C EO CARO 

100 WR ITE ( 6 » 600 3 N2»CT<LH> »KS»KT 
WRITE (6 f 610 3 i X ( 1 ) » ( XX ( J 3 » J=1 , 4 ) 

610 FORMAT ( 1H+ ► 23X > 14 » 25X » 4A6 3 
IFdXd) aGT . 0 ) GO TO 101 
LN-12 

GO TO 1000 

101 IF ( IX ( 1 ) .GT.10) Call EPSU) 

LN— 1 1 

GO TO 1000 

c E<3 CARD 

110 WRITE (6 » 600 ) N2»CTlLR) rKStKT 

WRITE (6 >611 3 (IX(j)»j=l#2)#(XX(J)fJ=l»7) 

611 FORMAT ( 1H+ t 23X* 214 # lbX»7F10. 3) 

LN~1 0 

60 TO 1000 

C OD CARD 

120 WRITE (6? 600 ) N2 t CT C LK 3 r KS , KT 
WRITE ( 6 o 612 ) IX 1 1 )'» C X X t U ) » J=1 ,4) 

612 FORMAT C 1H+ # 23X r 14 , 25X 1 4A6 3 
IFCIXC1) «GT « 0 ) GO TO 121 
LN=14 

GO TO 100Q 

121 I-F(IX(I) .GT.20) CALL EPSCl) 


LN=I3 

(30 TO 1000 

00 CARD 

130 WRITE (6 r 600) N2»CT(LR) *KS*KT 
WRITE (6*613) (XX< J) * J=l*3) 

613 FORMAT (1H+»23X* 314) 

UN™ 1 2 

GO TO 1000 

PD CARD 

140 WRITE {6*600) N2 * CT < lR ) *KS* KT 
WRITE ( 6 * 614 ) IX ( i ) # (XX(J) *J=1»6> 

614 FORMAT C 1H+ * 23X *14* 2bX * 6A6 ) 

IF ( IX ( 1) o GT « 0 ) GG To 141 

LN=16 

GO TO 1000 

141 IF(IXli) oGT.120) CALL EPS ( 1 ) 

IP=IX(1> 

LN"15 

GO TO 1000 

PQ CARD ” 

150 WR ITE ( 6 ► 600 ) N2 * CT (LR ) *KS * KT 
WR ITE ( 6 * 615 ) XX* l XX l j) * J=l*2) 

615 FORMAT ( 1H+ * 23X * 6 14 t 2F 10 . 3 ) 

00 152 J~1 * 4 

IF ( IX ( U ) .GT.50Q) 60 TO 151 

IF( (IX(J) *GT.O) .AND. (IX (U> .LE.20) ) GO TO 151 
IA0< IP# J)=IX(J) 

GO TO 152 

151 CALL EPSll) 

152 CONTINUE 

IF(XXCS) .GT.10) CALL EPS ( 1 ) 

IF ( IX (6 ) s GT • 500 ) CALL EPSC1) 

IF ( IX ( 6 ) ®LE ® 0 ) CALL EPS (3) 

LN=14 

60 TO 1000 

DD CARD 

160 WRITE <6 t 600) N2*CT <LR> *KS*KT 

WRITE (6*616) (IXIJ) * J~ 1 * 2 ) * ( X X ( U ) *0-1*4) 

616 FORMAT ( 1H+ * 23X * 214 * 21 X * 4 A6 5 
JD=IX(1) 

M= I X < 2 > 

IF ( IX ( 1 ) oGT«0) GO TO 161 
LN-2 

60 TO 10U0 

161 IF ( IX ( 1 ) «GT» 500 ) CALL EPS(l) 

IF( (M.GT.O) •AND«(M.LE.11> ) 60 TO 162 
CALL EPS ( 1 ) 

162 IF(ND.LT.IX(1)> ND=IX(1> 

LN-17 

GO TO 1000 

DQ CARD 

170 WRITE (6 *600 ) N 2 » C T ( L R ) * K 5 * K T 
WRITE (6*617) (XX (J) * J=1*M) 

617 FORMAT (8X*F10«2» 10F10»6) 


If C ABS ( SUM"! ® ) eSTeeOOOOOOl) CALL EPSUU 

172 MP=M+11 

WRITE ( 6 9 637 ) C XX ( J ) » j=12 » MP ) 

637 FORMAT (8Xt»iiElQ.2) 

LN-16 

IF CMaLQ « i ) GO TO 1001 
DO 173 J=2#M 

IF ( XX ( J+10 ) «LE . XX ( J+ll ) ) GO TO 173 
CALL LPS(6) 

GO TO 1001 

173 CONTINUE 
60 TO 1001 

c 2D card 

180 WRITE (6e6QG ) N2 » CT ( LR ) » KS > KT 

WR ITE (6 f 610 ) IX C 1) »(XX< J) » J=l#4) 

610 FORMAT (lH+»23X»I4»2bX» 4A6) 

IFUX(i) ,GT.O) GO TO 181 

LN=1 

60 TO 1000 

181 IF ( IX ( 1 ) .6T.120) CALL EPS C i ) 

LN=19 

GO TO 1000 

C 2C CARD 

190 WR ITE ( 6 # 600 ) N2»CT<LR) »KS»KT 
WRITE (6 » 619) IXU) » l XX < J) » J=1 » 4 ) 

619 FORMAT ( 1H+ ?23X»I4»20X»4F10«3) 

IF ( IX < I ) • 6T . 0 ) 60 To 191 
LN—18 

GO TO 1000 

191 IF(IXU) .6T.120) CALL EPS(l) 

LN=19 

60 TO 1000 

200 WRITE (6r 600 ) N2 > CT ( LR ) * KS * KT » N 

WRITE (6 t 605) ( IX ( J ) > J= 1 t 2 ) » ( XX ( J ) » J=1 » 4 ) 

f P — T X ( 1 ) 

IF(IP.GT.O) GO TO 201 
LN=4 

60 TO 1000 

201 IF ( IP. 61. 120) GO TO 2C8 
LK=IX(2> 

IFC (LK.LE.O) .OR* (LK.GT.500) > GO TO 58 
DO 202 J=l#4 

IF(XX(d)*XR(LK) .GT.l.) CALL FIPS (6) 

202 CONTINUE 
GO TO 209 

208 CALL EFS(l) 

209 LN=2U 

1000 WRITE (9) LR» IX»(XX(J) »J=1»7) 

IF (LRtNEoLX ) CALL EPS(7) 

RETURN 

1001 WR ITE ( 9 ) LR t XX 
IFCLR.NE.LX) CALL EPS(7) 

RETURN 

END 



SUBROUTINE EPS CEE) 

ERROR PRINTING SUBROUTINE 

COMMON/ T /LN » LR » I X ( 6 ) * * XX(22) , M » Nl * N2 * I ( 6 , 20 > *XC22*2Q) oKC22#20) * 

• KlL»NErKE(ll) * IAB<120»6> t DR ( 50 0 ) » XR ( 500 ) r NO t X0*AR(4) *KKT(20) *NP 

•LH (20 ) *KXT 
IF (LE ®EQ • 0 ) 60 TO 100 
NE-”NE^'I 

KE ( LE ) =KE (LL ) +1 

IF ( (KIL.EQ.l) .AND. CLE«.6T,7) ) RETURN 
NP=NP+3 

IF (NPoGE.55) CALL PTS 
WRITE (6 *600) NE 

600 FORMAT (/12H ******ERROR > 14 ) 

60 10(1*2*3*4*5*6* ?*8*9*10*li) * LE 
1 WRITE (6*601) 


601 FORMAT ( 1H+ * 17X * 30H INDEX OUT OF RANGE /) 

RETURN 

2 WRITE (6 *602) 

602 FORMAT (1H+* 17X»30hPREREQUISITF MISSING /) 

RETURN 

3 WRITE (6 #603} 

603 FORMAT ( 1H+ * X7X * 30HNECESSARY DATA MISSING /) 

RETURN 

4 WRITE (6 *604) 

604 FORMAT <lH+» 17X r 30HTOO MANY CARDS IN SET /) 

RETURN 

5 WRITE (6 *605) 

605 FORMAT ( 1H+ » 17X * 3 OH INDEX USED PREVIOUSLY /) 

RETURN 

6 WRITE (6 *606) 

606 FORMAT (1H+»17X#30hDATA OUT OF RANGE /) 

RETURN 

7 WR ITt (6 ? 60? ) 

607 FORMAT ( IH+» 17X» 30HCARD OUT OF SEQUENCE " /) 

RETURN 

Q WRITE (6*608) 

608 FORM AT(1 H+»17Xf3UH AFFECTED PART AREA ZERO /) 

RETURN 

9 WRITE (6 *609) 


609 FORMAT (1H+»17X»30HCONTACT AREA GT SURFACE AREA /) 
RETURN 

.10 WRITE (6 *610) 

610 FORMAT (1H+»17X»30HREMAINING AREA NEGATIVE ' /) 

RETURN 

11 WRITE (6*611) 

611 FORMA HlH+»m#3QHTOTAL PROB* NOT EQUAL To 1 /) 

RETURN 

100 WRITE (6*620) 

620 FORMAT <lHi//////2OX»27HMICRO0 I AL BURDEN DATA CHECK/) 
lF(NE.GT.O) GO TO 150 

WRITE (6 *621) 

621 FORMAT (24X*18HN0 ERRORS THIS RUN) 

RETURN 

150 W R 1 T E ( 6 * 6 2 2 ) 

622 FORMAT (20X»27H****'* ERROR SUMMARY ***♦*/) 

XOUT— 0 


• Hit s ; ■ <• A * * AJU * iU 

(XTE (6*623) KECu) 

)RMAT C 16 f 9H ERRORS -) 



156 WRITE(6»624) K£(U) 

624 FORMAT ( 16, 9H ERROR -) 

159 IOUT=l 

GO TO (101, 102, 103, 104, 105, 106, 107} , J 

101 WRITE (6,60 1 ) 

GO TO 152 

102 WR ITE ( 6 , 602 ) 

GO TO 152 

103 WRITE (6,603) 

GO TO 152 

104 WRITE ( 6 ,604 ) 

GO TO 152 

105 WR XTE ( 6 , 6Q5 ) 

GO TO 152 

106 WRITE (6,606 ) 

GO TO 152 

107 WR ITE ( 6 , 60? ) 

152 CONTINUE 
XF(KlL.NE.l) GO TO 154 
IF(KXT.EG.i) return 

IF ( I OUT oEQ.li CALL EXIT 
RETURN 

154 UO 153 J=8 » 1 1 

IF ( KE ( U > “1 ) 153,158,157 

157 WRITE ( 6, 623 ) KE(U) 

GO TO 160 

158 WRITE ( 6 , 624 ) Kt'(J) 

160 JM-J-7 

GO T0(108,109,110,lll) ,JM 

108 WRITE (6,608) 

GO TO 153 

109 WRITE(6,609) 

GO TO 153 

UO WRITE (6,610 ) 

GO TO 153 
ill WRITE (6,611 ) 

153 CONTINUE 
IF(KXT.EQ.l) RETURN 
IF(IOUT.EQ.i) CALL EXIT 
RETURN 

END 


subroutine pis 

PAGE TITLING SUBROUTINE 

COMMON/7 /LN»LR» IX (6) , XX ( 22 3 » Mf lv|i # N2 » I ( 6 » 20 ) » X ( 22 » 20 ) » K ( 22 r SO ) » 
.KIL»NE»Kt(U) »IAB(120»6) ?DR(500> » XR ( 500 5 t NO # 10 » AR ( 4 ) ? KKT<2Q) »NP» 

9 LH (20 ) f KXT 
WR ITE ( 6 » 600 ) 

600 FORMAT ( 115H1RECCRD t TYPE » S/TSK/ST 1 2 3 4 5 6 

.7 a 9 10 11 12 13/) 

NP~0 

RETURN 

END 



000 0 0000 


A I 6 

SUBROUTINE HPS 

c burden prediction subroutine 

C THIS PROGRAM READS FROM TAPE 9 ' 

COMMON KK* IK#KK»RUN(7) »KS'STG(7> 'KT'T$K(7) »DSC<6> » ID» 

.NE#AES(4) » lEC(lO) , IETQO) .AEC<10) ,AET<10) ,AED(10) ,AEF(l0f4) » 
«NO»IGT(20) » XOQ ( 20 ) » IQC (20 ) p 

• NP * DAC ( 120 f 2 ) ' IU ( 120 ) 'XAB(12Q»6) t AAG ( 120 ) » A AS (120 ) ? AAl ( 120 > » 

• NO rNX ( 572 ) * DR ( 572 Eli ) »XR(572»lD »K0# 

• NSfKNI (20) »IR(5) » AT? ( 4 ) » APA ( 2 ) . 

.NTrjT(lOO) » XMT(IOO) »XVTn.OO} »L1 »L2»L3»L4 »KTS» ITE 
DIMENSION XX ( 7 ) » IX(fa) 

3 FORMAT (34H1MICRCBIAL BURDEN PREDICTION MGDEL/4H RUN»I5»2H» »7A6) 
REWIND 9 
REWIND 11 

nt=o 

DO lb N=1 t 572 

15 NX(N)=-1 
DO 16 N-l ? 120 
DO 16 U=i»6 

16 IAB(N» J)=0 

100 READ (9) L»KK»IK»I1»I2»I5»I6»XX 
KK-1 INDICATES A NEW RUN? STAGE » AND TASK 
KK-2 INDICATES A NEW STAGE AND TaSK 
KK=3 INDICATES A NEW TASK 

KK»4 CALLS FOR A BURDEN DIFFERENCE DETERMINATION 
KK-5 CALLS FOR ZONE BURDEN WRlTE'OUT 

iK-o calls for zone definition data inputs 
ik-1 indicates zones correspond to parts 

KK=6 CALLS FOR A RESTART AND REQUIRES A BINARY DECK INPUT 
IF (KK »NL • 6 ) GO TO 102 

REWIND 2 A 

READ (2 ) KR t (RUN< J> ' J=l»7 ) »KS» (STO(J) i J= l» 7> » NF » ( AES ( J ) t J=l»4 ) , A 
. ( IEC ( J ) » IET ( J ) » AEC ( J ) , AET( J) » AED C J ) » (AEF( JcK) » K = 1 1 4 ) » J= 1 , NE ) » 
a NO' ( IOT (o> » 10Q ( J) » IOC { J) # J=1»N0) »NP» ( (DAC ( J e K ) t K=1 » 2 ) t IU ( J ) * 

• UAfcMJ'K) »K=1»6) »AAG(J> * AAS ( J ) ' AAT < J> ► J=1 » NP ) » ND » KO * NT t 
. ( JT ( J) * XMT ( J ) # XVT< J) # J=1»NT) »KTS» ID» ITE»Kl»K2 

101 READ (2) N»NX<N> # (OR(N»J> #XR(N»J> »J=l»ll) A 

IF (N »LT . 572 ) GO TO 101 
GO TO 100 

102 IF ( KK-3 ) 103 o 130 f 109 

103 IF (NT»EQ . 0 ) 60 TO iGo 

C 

C WRITE STACsE SUMMARY FOR PRECEEDlNG STAGE 

WR ITE ( 6 » 3 ) KR » (RUN( J) » J=l»7) 

WRITE (6»4)KS» <STG( J) » J=l» 7) 

4 FORMAT ( bH STAGE' I3»2H» '7A6/15H STAGE SUMMARY-// 

• 3X»4HTASK#7X»4HMEAN»6X» 6HF INISH/ 13X » 6HBURDEN * 6X» 4HT IME/ ) 

DO 104 1 = 1 » NT 

104 WRITE (6» 5) JT ( I ) » XmT ( I ) » XVT ( I ) 

5 FORMAT (X7»E13.3»F10*3) 

106 IF ( KK « GT » 0 ) GO TO 109 
IF (K2 • EU • 3 ) GO TO 108 

WRITE (2) KR' <RUN( J> r J=l»7> »KS* (STG< J> * J=l»7> »NE» (AES< J) $ U 1 '■* * r •’ ? 

. (IEC(J) 'IET(u) 'AEC(J) s> AET ( J ) * AED ( J ) » ( AEF ( J 'K ) ' K=1 » 4 ) » J= 1 * NE 3 t 

• NO? ( IOT ( J ) ' IOQ ( J ) ' IOC ( J ) t J=1»N0> ? NP > ( (DAC (u'K) »K=1'2) f IU( J) » 

• ( IAB(J'K) »K=1»6) 'AAG(J) ' AAS ( J ) 'AAT(J) r J=1 » NP ) » ND t KC » NT » 

• (U* v *JI ' XMT ( U ) t XV i » -J ) > J= 1 ' NI ) f KTS» ID' ITE»Kl*K2 

N— 1 

107 IF ( NX ( N ) ® GE * 0 5 WRITEC2) Nr NX( N) ' (DR ( N t J ) # XR ( Nr J) »J=1 r 1 1 ) 

N=N+1 

IF (N *NE s 572 } GO TO 107 


c <t «3C 



o o n o o n o no 


A: 7 

WRITE. (2) N»NX(1) » (OR(h J) *XR(1»J) t J=l?ll) 

iOfl M~0 ' 

WRITE (11) M#KR»KS»M»M» (DSC(J) » J=l»4) (RUN ( J> » J=le7) » C StG C JS » 

» U“ 1 » 7 ) * { JJ { j } t J— 1 1 la ) 

RETURN 

X 09 GO TO <110»120»130»140*150) »KK 
C 

C READ RUN NUMBER AND DESCRIPTION 

110 READ (9) l.*KR» (IX(J) *J=1*5> »RUN 
WRITt (6*600) KR» (RUiM(J) » J=l*7) 

600 FORMAT (34HlMICR0bIAi_ DURDEN PREDICTION MQDEL//4H RUNpI2p2X^7a6//) 

ID=0 
MD = 0 
iU~;i 
1 

NE = 0 
N0=0 
NP=0 
ND"0 

ns=o 

NT~0 

DO 111 1=1*120 
IU(I)=0 
AAG ( I ) =0 9 
AAS ( I ) =0 o 
AAT ( I ) =0 * 

00 111 J=1 

111 1 AB ( I * U ) -0 

UO 112 N=1 t 572 

112 NX ( N ) -”1 
C 

C READ STAGE NUMBER AND DESCRIPTION 

120 READ (9) L»KS*(XX( J) * J=l»5) *STG 
NT = 0 
M=1 

WRITE (11) M * KR * KS * M o M * ( DSC ( J ) * J=1 » 4 ) * M * M * (RUN(J) r J=l»7) * (STG( J) * 

. J=l» 7) * (uT ( J ) * J=l*la) 

READ TASK NUMBER AND DESCRIPTION 
130 READ (9) L»KT»L1»L2»l 3»L4»IX(1) » TSK 
WRITE, (6*612) KS*STG*KT*TSK 
612 FORMAT ( 6H1STA6E * 13 e 2X * 7A6/5H TASK * 14 ? 2X* 7 a6/ 5 
LI AND L2 ARE PREREQUISITE TIME DISTRIBUTIONS 
L3 IS THE FINISH TIME DISTRIBUTION IF NEEDED ’ 

AS PREREQUISITE FOR ANOTHER TASK 
L4 IS THE INDEX FOR SAVING THE TOTAL MICROBIAL BURDEN 
IF NEEDED FOR USE IN DETERMINING A BURDEN DIFFERENCE 

DETERMINE TASK START TIME 

KTS=0 

IF (LI *LE 9 0 ) GO TO 134 
LI=Ll+530 

IF (L2aG! ®Q)GQ TO 132 
KTS=L1 
GO TO 134 
132 L2=L24-530 

CALL HCS ( i, a t L2 p KTS’* 5 ) 

IF ( 10 ol-E • 0 ) GO TO 134 

IF C NX ( KTS ) o 61 9 NX C ID ) ) CALL HAS ( KTS ) 

IF C 1 1 o 6 Y o 0 5 K 1 = 11 


134 
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f \ ! b 

IF ( X2#6T ®q)K2=X2 
IF ( I5.LE • 0 ) GO TO 138 
IF( ID.GT.O) MO=NXUO) 

ID-15 

C RESTORE THE FIRST 20 DISTRIBUTIONS ' 

00 13? 1-1*20 
0— I 4552 

IF(nX(J) .LT.O) GO to 13? 

CALL HES (1*0) 

137 CONTINUE 

138 IF( Ib.GT.O) ITE=Ib 
NT=NT+1 

CALL MbS(M0*Kl*K2) 

1F(K2.NE.2) GO TO 100 
REWIND 2 

WRITE (2) RR * (RUN(0> *0=1*7) *KS* (STG( J) *0=1*7) * NE * (AES( J) *0=1*4) * 

. (1EC(0) *IET(0) *AEC(o> *AET( 0 ) »AED(0) * (AEF(0*K) *K = 1*4) *0=1*NE) * 

. NO* (IOT(u) * IOQ(O) *10C(0) *0=1* NO) »NP» ( (DAC(0*K) * K=1 * 2 ) * IU ( 0 ) # 

. ( I AB ( J * K ) *K=1 *6) * AAG(O) * AAS(j) >AAT( j) * 0=1 * NP ) * ND * KO * NT * 

. ( JT(O) * XMT l J ) * X V T ( 0 ) *0=1* NT) * KTS * ID * ITE * Kl * K2 
N=1 

139 IF (NX(N) • GE • 0 ) WRITE(2) N #NX ( N ) * (DR < N * 0 ) * XR t N * 0) * 0=1 * 11 ) 

N=N+i 

IF (N • Nt s 572 ) GO TO 139 

WRITE (2) NfNX(l) * (DR C 1*0) »XR(1*J) *0=1*11) 

GO TO 100 

DETERMINE AND WRITE BURDEN DIFFERENCE 

140 WRITE ( 6 * 3 ) KR * (RUN(O) *0=1*7) 

WR ITE ( 6 * 60 1 ) 

601 FORMAT ( //18H BURDEN DIFFERENCE/) 

L1=I 1+530 
L2=I24530 

WRITE (6*602 ) I1*DR<L1*1) 

602 FORMAT ( /17H BURDEN STORED IN* I3*7H* AREA= * F9 . 2 * 4X * 6HBURDEN ) 

CALL MWS(Ll) 

WRITE (6*602) I2*DR(L2*1) 

CALL HWS ( L2 ) 

A1=DR(L1* 1) 

A2=0R ( L2 * 1 ) 

DR(L1*1)=AVF(L1) 

DR(L2*1)=AVF(L2) 

CALL HCS(L2*Li*55l*2) 

CALL HAS ( 551 ) 

WRITE lt> * 601 ) 

CALL HWS (551 ) 

DR(L1* 1)=A1 
DR(L2*1)=A2 
60 TO 100 

DETERMINE AND WRITE ZONE BURDEN DISTRIBUTIONS 
150 WRITE ( 6 ► 3 ) KR * ( RUN ( 0) *0=1*7) 

WRITE (6*7)KS»KT 

7 FORMAT (43H MICROBIAL BURDEN BY ZONES (FOLLOWING STAGE*I3*bH* TASK* 
. 1 3 * 2 H ) " / ) 

IF ( IK ® EG) • 1 > GO TO 180 

USE ZONE DEFINITION INPUTS 
155 READ (9) L * XZ * (IX < 0 ) * 0=1 * 5 ) * DSC » XX ( 1 ) 

IF ( IZ • LE • 0 ) 60 TO 100 
A— 0 * 


CALL HES ( 510 » KO ) 

160 READ (9) L» IP» ( IXCJ) » J=l»5> » AR» <XX( J) r J=l»3) - 

c ip is the Part 

C AR(J) THE FRACTION OF SURFACE J OF PART IP BELONGING TO ZONE IZ 

tc / to i r ffi 'T'l 

jrvAr «Ul«u<v)U f U ' U 

DO 165 J=i?4 

IF ( AR( J) »lQ»Os)GO TO 165 
IFSIABdPf J) .LQ.UJGO TO 165 
F=AR(J) 

IB=IAB(1P»«J> 

CALL HMS ( XB e F e SOS ) 

A=A+DR ( 509 » 1 ) 

DR(509»1)=AVF(509)- 
CALL HCS(510*509»510»1) 

165 CONTINUE 
GO TO 160 

170 WRITE ( 6 » 9 ) IZ » (DSC ( J) » J=l»4) * A 

9 FORMAT <5H ZONE » 14 F2X » 4A6 » 7H t AREA= » F8 . 3/33X t 1 1HZONE BURDEN) 

CALL HWSC51U) 

GO TO 155 
C 

C USE PARTS AS ZONES' 

160 DO 190 I=1»NP 

IF ( IU( I) oEQ b 0 ) GO TO 190 
A--0 » 

CALL HES ( 510 * KO ) 

DO 185 0=1 t 4 
IB=IAB(I»o) 

IF(lB.LEeO) GO TO 165 
A=A+DR( IBrl) 

CALL HES ( 509 ? IB) 

DR(509*1)=AVF(IB) 

CALL HCS(510»509f 5lQ#l) 

165 CONTINUE 

WRIT^(6*10) I»CAC(I»1) » DAC ( I e 2 ) * OR (510*1) 

10 FORMAT (5H ZONE # 14 , 2X » 2A6 t 7H » AREA=»F7.2» BH» BURDEN) A 

CALL HWS ( 510 ) 

190 CONTINUE 
GO TO 100 
END 
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SUBROUTINE MBS ( MD . K1 . K2 ) 

C MICROBIAL BUILDUP SUBROUTINE 

COMMON KK.IK.KR.RUN17) #KS.ST6(7> .KT.TSM7) #DSC(6> »XD. 

•NE# AE8<4> » IEC (10) . IET(IO) » AEC ( 10) * AET ( 10 ) , AED(IO) .AEF ( 10.4 ' » 

•NO. IOT (20) » 100(20) » IOC (20) . 

•NP»QAC(12U»2) . ZU ( 120 ) » X A6 ( 120 . 6 ) .AAG(120) rAAS(120) »AAT( 120 > » 

• NU » NX ( 572 ) . DR ( 572 . 11 ) . XR 1 572 » 1 1 ) » KO . 

•NS»KKT(20) . IR(5> #AR(4> »APA(2> » 

• NT. JT(IOO) .XMTC100) »XVT(lQ0).Ll.L2»L3.L4.KTS»lTE 
DIMENSION iX(fo) » XX ( 7 ) 

GO TO (100.120.140.160.200) . IK 
C 

C ENVIRONMENTS INPUTS" 

100 READ (9) L » IX » AES » <XX< J) . J=l»3) 

WRITE (6. 8) ( AES ( J ) » J=l»4> 

8 FORMAT ( //30H ENVIRONMENTS INPUTS /5h AES=»4E12.4) 

C AES ( *J ) IS THE SURFACE LIFETIME MODIFIER FOR SURFACE J 

DO 110 1=1.10 

READ (9) L»N»(IX(J).J=1»5) »DSC»XX(1) 

IF ( N • LE » 0 ) GO TO 120 

READ (9) L»IEC(N> »IET(N) . <IX(J) .J=l»4> »AEC(N> .AET(N) »AED(N> » 

. (AEF(N» J) » J=1.4) 

C IEC » AEC DESCRIBE THE ENVIRONMENTS BIOTA CONCENTRATION 

C IETfAET DESCRIBE THE REFERENCE ACCRETION TIME 

C AED IS THE AIRBOURNE CONCENTRATION PER MAN 

C AEF IS THE RATE AT WHICH. FOR ALL OTHER FACTORS STANDARD# 

C 8I0TA REACH EACH SURFACE 

WRITE <6.20 >N» (DSC (J) »J=1* 4) . IEC(N) » AEC (N) . IET(N) » AET (N> .AEO(N> . 

. ( AEF ( N » J ) » J=1.4) 

20 FORMAT (/ 14 #1X»4A6,4X,4HIEC=# 13# 6H* AEC= » F7 . 4 , 6H . IET=. I3.6H. AET- » 
•F7.2.6H. AED=#F7.h#6H# AEF=#4F6.2) 

110 IF(NE.LT.N)Nt=N 

REAO (9) L»N» <IX(o> »J=1»5) »DSC»XX<1> 

C 

C OPERATIONS INPUTS- 

120 WRITE (6.9) 

9 FORMAT ( //3GH OPERATIONS INPUTS 
DO 130 1 = 1 » 20 

READ (9) L#N# ( IX( J) # J=l#5) #DSC#XX( 1) 

IF (N »LE • 0 ) GO TO 140 

READ (9) L#IOT(N) . I0Q<N> . IOC ( N > » ( IX ( J) . J=1 . 3 ) , XX 
WRITE ( 6 # 21 ) N # (DSC ( J) » J=1 » 4> » IOT <N> » IOQ(N> » IOC (N> 

21 FORMAT (/I4#lX»4A6r4X#4HI0T=» 1 3 » 6h » IO<2=» I3»6H» I0C=» 13) 

C IOT IS THE OPERATION TIME INTERVAL 

C IOQ IS THE DIRTINESS FACTOR 

C IOC IS THE BIOTA CONCENTRATION FOR CONTACT CONTAMINATION 

130 IF(NO.LT.N>NO=N 

READ (9) L»N»<IX(j)*J=l»5) .DSC.XX(l) 

C 

C PARTS INPUTS" 

140 WRITE (6# 10 ) 

10 FORMAT ( //30h PARTS INPUTS 
DO ISO 1=1.120 

READ (9) L » N . ( I X ( J ) » J=1 . 5 ) » DSC . XX ( 1 ) 

IF (N.LE.O)GO TO 160 
UAC ( N . 1 ) =DSC ( 5 ) 

UAC ( N » 2 ) =DSC ( u ) 

IU(N)=0 
AAT ( N ) =0 • 

IF(KTS.GT.O)AAT(N)=DR(KTS» 1) 

READ (9) L» CIABCN. J) » J=l»6) #AAG?N) .AAS(N) r (XX(J> . J=1.5) 
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IF < I Ab ( Ne 6 ) «GT t, 0 ) GO TO 142 
AAG(N)= 0 « 

AAS(N)=0« 

142 CONTINUE 

UAC XS THE PERMANENT ALPHAMERIC DESCRIPTION OF PART N 
1A8(J> INDICATES THE DISTRIBUTION FOR AREA/BURDEN FOR 
J=l# TOP SURFACE 
J~2 1 OTHER EXTERIOR SURFACE 
J-3 p MATED SURFACE 
U=4 o OCCLUDED SURFACE 

IAB(5) IS THE ENVIRONMENT INDEX (MAY BE LEFT BLANK) 

1 AB ( 6 ) IS THE RETENTION DISTRIBUTION FOR CONTACT 
AAG IS THE RETENTION DISTRIBUTION MEAN FOR FALLOUT 
A AS IS THE RETENTION DISTRIBUTION MEAN FOR CONTACT 
AAT IS THE LAST TIME OF ACCRETION UPDATE 

WRITE ( 6 » 22 ) N » (DSC(J) , J=1 1 6 ) > ( I AB ( N » J ) » J= 1 » 6 ) > AAG(N) »AAS(N) 

22 F0RMAT(/I4»lXr6A6»4X»4HlAB=*6I4*6H» AAG= » P? . 3 , 6 H t AAS=*F7 6 3) 
ISO IF(NP.LT.N)NP=N 

READ (9) L>N> ( IX ( J) t J=l*5) jQ5C»XX(l> 

DISTRIBUTIONS INPUTS - 
160 WRITE ( 6 » 11 ) 

11 F0RMATC//30H DISTRIBUTIONS INPUTS -//) 

00 17S l=ie5G0 

READ (9) L»N»M» (IX(J) »J=1»4) »nSC»XX(l) 

M=1 indicates the distribution is a constant 

IF ( N • LE • 0 ) GO TO 179 

READ (9) l» <DR(N* J) » J= l» 11) * (XRCN* J) »J=l»ll) 

NX(N)=M 

IF (OR (N» 1 ) . LE » 0 a ) DR(N»1)=AVF (N) 

IFCN.GToSO) GO TO 171 
«j=N+552 
CALL HES(JpN) 

171 IF ( ID.IE.G) GO TO 172 
if ( m • gt • nx ( xo ) ) call haS(n) 

172 WRlTt(6»23)N» (DSC(J) »J=1»4) »DR(N»1) 

23 FORMAT! X4»1X#4A6/5X» 12HC0EFFICIE'NT=»E11.4»6X» 12HDISTRIBUTI0N) 
CALL hWS ( N ) 

XF(ND.LT.N) nd=n 

175 CONTINUE 

READ (9) L,»N*M»'<IX(J) »J=1»4) »CSC»XX(1) 

179 IF( ID.LE.O) GO TO 190 

IF (MO » E0-) NX ( ID) ) GO TO 190 
DO ISO I=l»20 

IF C NX < I ) .GT.NX(ID)) CALL HAS ( I ) 

160 CONTINUE 

IF ( MD ® LE • NX ( ID ) ) GO TO 189 
DO 185 1=21 o 552 

IF ( NX ( I ) 9 GT 9 NX ( ID ) ) CALL HAS ( I ) 

165 CONTINUE 
169 MD-NX(ID) 

THE ZERO DISTRIBUTION 
190 IF (NO. GE. 500) WRXTE(6»630) 

630 FORMAT C///33H **** TOO MANY DISTRIBUTIONS ****///) 

ND=ND+1 

KO=ND 
NX CKO ) =1 

DR(KO»l)=0 o 
X R C K 0 f 1 ) = 0 o 

c subtask level microbial buildup determination 
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200 IF(KTS»EQ®0)KTS-K0 
€ KT5 IS THE TASK START TIME 
LtsKTS 
NS=0 

00 295 1-1 s 20 

READ (9) L»N*N1»N2» UX(U) * J=l»3> *0SC*XX(1> 

C Nl AND N2 ARE PREREQUISITE 5URTASKS 
IF ( N . LE « 0 > GO TC 300 
WRITE i6» 12) N* (DSC ( J) *J=I*4> 

12 FORMAT ( /8H SUBTASK* I3*2H* »4A6> ' 

IF (NSsLTeN) NS-N 
IF (Nl »GT . 0 ) GO TO 205 
I CT-KTS 
GO TO 215 
205 J1=KKT(N1) 

IF (N2.GT.O)GO TO 210 
ITT=J1 
GO TO 215 
210 J2=KKT(N2) 

ITT-5 10 

CALL HCS ( Jl * J2 * ITT » S ) 

C ITT IS THE SUBTASK START TIME INDEX 

215 T1=UR(ITT*1> 

M=2 

WRITE (ID M*KR*KS*KT*N* (DSC( J) * J=l»4) * ITT * NX ( ITT) * ( DR ( ITT * U ) * J=l# 
.11) * (XR< ITT* J) »U=1* 11) * <AEC( J) *J=1*5) *N1*N2* ( JT( J) * J=l*3> 

NC = 0 

C ENVIRONMENT* AREA/BIOTA BURDEN CHANGES- 

220 READ (9) L*K* IR* AR* (XX< J) * J=l,3) 

K IS THE CHANGE IDENTIFIER 

I R ( 1 ) * I R ( 2 > ARE THE NEW PART AND SURFACE INDICES 
IR ( 3 ) * IR ( 4 ) ARE THE CONTRIBUTING PART AND SURFACE 
IR ( 5 ) IS THE NEW ENVIRONMENT (IF ANY) FOR PART IRQ) 

AH ( 1 ) IS THE AREA CHANGED TO IR(1)*IR(2) FROM IR ( 3 ) * IR(4) 

AR ( 2 ) IS THE NEW AAG FACTOR 
ARC 3) IS THE NEW AAS FACTOR 

IF IR ( 5 ) * AR ( 1 ) * AR ( 2 > * AR ( 3 ) ZERO* NO CHANGE IS MADE 
IF ( K *LE • 0 ) GO TO 240 
NC=NC+1 

IF(NC.EQ.l) WRITE(6*16) 

16 FORMAT ( /toX * 42HENV IRONMENT/ ARE A/RETENT ION FACTOR CHANGES-/) 

1 1 = I R ( 1 ) 

IFCNP.LT. I1)NP=I1 
IF ( AR ( 1 ) .LE.0.)G0 TO 232 
CALL MAStllfTl) 

IU(I1)=IU(I1)+1 
I2-XR ( 2 ) 

13=1R(3) 

CALL MAS ( 13 * T1 ) 

IU(I3)=IU(I3)+1 
I4=XK(4> 

1B=IAB( II f 12) 

IA=IAB(I3»I4) 

IF ( 16 .GT . 0 ) GO TO 2'22 
IFCND.6E.50Q) WRITE(6*630) 

ND=ND+1 
IB=ND 

IAB ( i 1 * 12 ) = XB 
CALL HES ( IB *KO ) 

222 IC=NU+i 

IF CIRCS) I A8 ( II » 5) -IR ( b ) 
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IF ( AR (2 > »©T * 0 s } AAG ( 1 1 ) ~AR (2 ) 
IF(AR(3) .GT.O.) AAS(Ii)=AR(3) 
IF(DR( IA» 1) .LE»0. )'GQ TO 220 

F=AR(1)/0R( IA» i) 

IF(F«LT« .995 ) GO TO 230 

f=i. 

IAB( 13* I4)=0 
230 CALL HMS(IA»F» IC) 

A=DR(113»1)+0RUC»1) 

DR(IB»1)=AVF(IB) 

DR(IC»1)=AVF(IC) 

CALL HCS( XB* lCr IB» 1) 

0RQb»i)=A 
F= l.-F 

CALL HMS( IA»F» I A > 

GO TO 220 

232 IF ( IR ( 5 ) » GT . 0 ) IAB(Il*b)“XR(5) 
IF< AR(2) .GT.O. ) AAG C 1 1 ) "AR < 2 ) 
IF ( AR ( 3 ) » GT o 0 9 ) AAS( I1)=AR(3) 
60 TO 220 


OPERATIONAL LEVEL BUILDUP- 
240 1T=5!0+N 

CALL HLS ( IT » ITT ) 

ITT IS THE SUDTASK START TIME DISTRIBUTION INDEX 
1T=510+N IS THE SUBTASK CURRENT TIME DISTRIBUTION INDEX 
KKT(N)=IT 
I E“ ITE 

250 READ (9) L * 10 # IKE * ( IX ( J) » J= 1 » 4 ) *AKT * AKQ * ( XX ( J ) » J=1 » 5 ) 

10 IS THE OPERATION 

IKE IS THE OPERATION ENVIRONMENT WHEN DIFFERENT FROM 

THE TASK ENVIRONMENT OR THE PREVIOUS OPERATION ENVIRONMENT 

akt is the operation time modifier 

AKQ IS THE NUMBER OF MEN 

10 LT 0 INDICATES A DECONTAMINATION OPERATION FOR WHICH 
AKT IS THE FIXED OPERATION INTERVAL 
Xf(XQ)252 *290*258 


DECONTAMINATION 


252 IO=-IO 

T1=DH( IT» 1) 

IF ( AKT «LE • 0 » ) 60 TO 253 

X7"ND^7 

NX (17) si 

DR(I7#1)=AKT 

XR(I'7*1)=AKT 

CALL HCS ( IT * X7 * IT * 1 ) 

253 XF(K1.GE.2) WR XTE ( 6 » 37 ) 10 * Tl * OR < IT * 1 ? 

3? FORMAT <5X»9H0PERAT ION r 13 »25 h» (DECONTAMINATION)? FROM*F?*2*3h TO? 
• F? s2?6H HOURS) 

T1=DR(IT»1) 

254 READ (9) L » XP »LK » ( IX ( J > t J=1 » 4 ) t AR » ( XX ( J ) * jsl t 3) 

IP IS THE PART AFFECTED 

LK IS THE CURVE DESCRIBING FRACTION Of BIOTA REMOVED 

AR( J) IS THE MEAN FRACTION OF BIOTA REMOVED FROM SURFACE J 

XP ( XP * LE » 0 ) GO TO 250 

IU(IP) = IUMP)41 ••• • • • * ■ .... 


QT=T1-AAT(IP) 

CALL MAS ( XP * Tl ) 

DO 256 J=l»4 

*s- e** n * t * % * *e /** f\ gt 

i r \ A n > <J I e LC. ® 41 o / \3 U i w C DO 


A 


i B~ 1 AB < XP#«J) 

IF( 16 .LEe 0) GO TO 256 
A~DR ( IB # 1 ) 

IF(A.LE.O» > GO TO 256 
IF CDT®LE»Oe ) GO TO 255 
M— 8 

WRITE (ID M#KR#KS#KT#N# (DSC ( L > # L=1 # 4 ) #X8#NX<IR> » ( OR ( IB # L > #L»1 # 11 ) 
. ( XR ( 10 # L ) #L=1»11) #A#T1» { DAC ( IP » L) # L=1 » 2 ) # A # IP # J # 10 » I AB < IF » 5 > , IT 
255 F=(1.-AR( J) )/DR(LK'» L) 

IC=ND+1 

CALL HMS<LK#F»IC.) ' 

DR(IB»1)=AVF(IB> 

UR(IC#i)=AVF(IC) 

CALL HCS(IB» IC# iBr 3) 

BQ=DR(IB»1) 

DR( IB# U=A 
M=9 

i«R ITE ( 1 1 ) M» KR»KS#KT»N» (DSC(L) #L=1#4> » IB # NX ( IB > » ( DR ( IB # L ) t L=s 1 , 1 1 ) 
• (XR(IB»L) »L = 1» ID »A»T1» AKT#AR( J) » A » IP r J » 10 #LK » IP 
1FSK1„£G«3> WRITE (6, 39) IP, DAC ( IPr 1> #DAC( IP# 2 ) » J » A # AR ( J ) t BD 
39 FGRMAT(10X#4HPART#I4»2X»2A6»9H# SURFACE# I2»7H» AREAs#F8.3# 
e 19H # FRACTION REMOVED^ »F5. 3 # 17X # 9H » BURDEN- » £10 . 3 ) 

256 CONTINUE 
GO TO 25H 
C 

C FALLOUT CONTAMINATION- 

258 lF(IKL®GTaO)IE-IKE 
I2-N0+2 
I3=NLH3 
I4=NQ-^4 
I5-ND+5 
I6=NU+6 
I7SN04-7 

ic=iocdo) 

10=100(10) 

II=IET( IE) 

JC=IEC(IE) 

ITK=I0T(IG> 

T 1=QR ( IT # 1 ) 

IF(ITK.EO.O)GO TO 259 
F = AKT/DRUTK»1) 

CALL HMS ( ITK # F # 17) 

0R(I7#I)=AVF(I7) 

CALL HCS( IT# 17# IT#1) 

C 17 INDICATES THE OPERATION TIME DISTRIBUTION 

259 IF(K1.GE.2> WRITE<6#29> 10 » IE # T1 # DR ( IT » 1 ) » AKQ 

29 FORMAT ( 5X » 9H0PER AT ION# 13# 13H# ENVIRONMENT # 13#6H# FR0M»F7.2»3h TO# 
« F7 , 2 # 7H HOURS# #F3.0»4H MEN ) 

AQ=AKQ/QR(XG#1) 

Q x AQ*AED ( IE) 

CALL HMS ( IG# G# IS ) 

0R(I5»1)=AVF<I5> 

C 15 IS THE CURVE G*t) 

AC=AEC(IE)/DR(JC»1 ) 

CALL HMS (JC# AC# 16) 

DRC I6»1)=AVK( 16) 

C 16 IS THE CURVE C 

CALL HC S i 1 b # 1 6 f 1 6 1 .1 ) 

c ie indicates the total fallout source concentration (c+o«d> 

260 READ (9) L # IP » LS » ( IX ( J) # J=1 » 4 5 # APD # APC # APS # APA # ( XX ( J ) » J=1 » 2 ) 

C IP IS the part affected 
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L5 IS THE TOOL STICKINESS DISTRIBUTION 

apd is not used in this version 

arc is the tool biota CONCENTRATION modifier 

APS IS the mean tool stickiness ' 

apa is the contact area for each touched surface 

IF ( IP »LE « 0 ) GO TO 250 
IU(IP)SIU(IP)+1 

UTS j 1 “hA ! i IP i 

CALL MAS ( IP » T1 ) 

IF ( AO . LE # 0 « ) GO TO 275 
IF ( It • LE • 0 ) GO TO 275 
IF ( ItK »LE • 0 ) GO TO 275 
AAT(IP>=DR(IT»1> 

DO 270 J=1 »4 
IB = IAB ( IP p J ) 

F=AEF(IE»J)*AAG(IP) 

IF ( IB .LE o G ) GO TO 270 
IF(F.LE.O* )GO TO 270 
A=DR<IBf X) 

IFCA«lE 6 O 0 )GO TO 270 
IF CDT « Lt « 0 a ) GO TO 2<a3 


C 


WRITE (11) M ( KR » KS t KT t N t (DSC(L) »L=1»4) * IB»NX( IB> » <DR(XB»L) »L=Irll) t 
. (XR(Xd’L) rL=lrll) »A»Tir (DAC(IP#L) » L-I » 2 ) r A * I P » J HO t I AB UP » 5 > e IT 

263 .F=F*A 

CALL HMSU6»F»I5) 

15 IS THE CURVE A*R = A*F*G* ( C+Q*D ) 


C 

C 

c 


265 


C 

c 


270 


EE=AET (I E > * AES ( J ) / DR ( 1 1 » X ) 

CALL HMS ( 1 1 p EE » 12 ) 

UR(I2»1>=AVF(I2) 

12 IS THE CURVE V 
CALL HCS(l5fI2»l4,3) 

14 |S THE CURVE A$V*R 
CALL HCS ( I 7 t 12 » 12 t 4 ) 

12 IS THE CURVE T/V 
M=NX(I2) 

NX( I3)=M 

DO 265 JJ=i»M 

DRC I3» JJ)=OR( I2» JJ) 

XR<I2» JJ>=EXP(-XR( I2» JJ) ) 

XR(I3» JJ)=l.-XR(I2»Jo> 

DR(I2»1)=AVF(I2) 

URU3rl)=AVFU3) 

12 IS THE CURVE EXP < -T/V) 

X3 xs the curve x-exp(-t/v> 

CALL HCS(I4» I3»14»3) 

OR(IB» 1)=AVF(IB) 

CALL HCS(I8»I2»I2#3) 

CALL HCS<I2» 14 o IB»1) 

CALL HAS (IB) 

BD-DH ( IB » X ) 

DRC IB»1)=A 
M=3 

WRITE ( XX ) MeKR?K5»KT*N# (DSC(L) »L=1»4) rlBrNX(Iff) » (DR( IB»L> »L=X»U) » 
( XR ( IB »L ) #L=X» ID f A» Tl»AKT»AKQ*Ar IPr J> IO»lEr IP 


IF r (K1 oEG a 3) WRITE (6 » 30 ) XP#OACC IP#1> t DAC ( IP * 2 ) t J » A r APC » APA ( J ) ?8D 


FORMAT ( J. OX 
14 m. ? TOOL 

CO NT 1 NUE 


4MO&P 


w f\ i 


14 f 2X » 2A6 » 9H t SURFACE »12»7H» AREA= »F8.3» 
Fo«0»15H» AREA TOUCHED 2 » F6 * 3 * 9H t BURDENS# EiO . 3 ) 


L 

C 


C 0 N T A C I C 0 N T A M I M A T 1 0 N - 
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275 IF (IC.IE.Q)GO TO 2'60 
DR(IC»1)=AVF(IC) 

IG=IAB< IPrb) 

IF(IGoLE.O) GO TO ‘260 A 

DR( IGr 1 ) =AVF ( 1G > 

Sl= e 5*AASC IP)*APC 
IF (LS»LE *0)60 TO 2.60 
DR(LS» 1)=AVF(LS) 

IF(OR(LS» i> *IE 8 U 6 )G0 TO 260 
S2=.5*APS/DR(LS'l) 

B-S1*S2 

IF(B«EQ.O. )GC TO 260 
S1=S1/(DR(XC'1)*DR(IG»1) ) 

DO 260 J=l# 2 
IB-I AB ( IP i J ) 

IF ( IB • LE • Q ) GO TO 260 
A-GR ( IB ' 1 ) 

DRUB' 1)=AVF (10) 

IF ( A ®LL » 0 e ) GQ TO 280 

IF ( APA ( U ) aLEeOa ) GO TO 280 

F-Si^APA ( J ) 

CALL HMS(IC»F#I4) 

CALL HC5( I4» IG* I4#3) 

F=S2*APA(J)/A 
CALL HMS ( LS * F » 15 ) 

CALL HCS(ia#I5»I5»3) 

CALL HCS ( IB ' 15 » IB » 2 ) 

CALL HCS(IB#I4»IB»1) 

CALL HAS (IB) 

OR ( IB c 1 ) = A 
M=4 

WRITE ( H ) M'KR'KS'KT'N' (OSC(L) »L=1»4) ' IB ' NX ( IB ) > (DR ( IB ' L ) > L=1 ' 11 ) ► 

. (XR(IB'L) »L = 1» U> » A t APA ( J ) 'AAS(IP) t APS' APC » IP » J' I A8 ( IP ' 6 ) »IC»IP 
260 CONTINUE 
GO TO 260 

290 IF(DR(LL'1> .LT.OR(ITd) )LL=IT 
295 CONTINUE 

READ (9) L»N»N1»N2' ( IX( J) » J=1'3) 'DSC'XX(l) 

300 IF (L3 8 LE b 0 ) GO TO 302 
L3=L3+53U 
CALL HES(LSiLL) 

302 CONTINUE 
T=DR (LL* 1 ) 

WRITE (6»60l ) KR ' KS » KT ' (TSK(J) » J=1'7) 

601 FORMAT (13H1TASK SUMMARY/4H RUN'I2'7H' STAGE' I3'6H» TASK » 14 » 2X » 7A6 » 
®//29X 1 15HTASK START TIME) 

CALL HWS(KTS) 

WRITE ( 6 ' 604 ) 

604 FORMAT (29X» l&HTASK FINISH TIME) 

CALL HWS(LL) 

WR ITE ( 6 ' 60S ) 

b05 FORMAT (/36H BURDEN DISTRIBUTIONS AT END OF TASK//50H BURDENS BY ZO 
• NE and SURFACE 
DO 303 J-1'4 
UJsSOO^u 
CALL HES(JJ'KO) 

303 AR(J>=0. 

00 310 1=1 'NP 

IF ( IU ( I ) ,F_Q.O)GO TO 310 

CALL MAS(I'T) 

WRITE (6 '602) I# (DAC(I'J) » J=l»2) 
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602 FORMAT (/SH ZONE t 14 * 2X # 2A6) 

DO 305 

JJ=5Q0*J 
IB=IA}J(Ip J) 

IF ( ia.LE.0>&0 TO 305 

IF (DR ( ID* 1 ) • LE • 0 » ) GO TO 305 

A“DR ( IB * 1 ) 

AR( J)=AR( J)+A 
0R(IB»1)=AVP(IB) 

CALL HCS ( OvJ * IB > UU > 1 ) 

DRC IB* D=A 

WRITE (6*603) J»QR(lD*l) 

603 FORMAT (5X» l4HSURFACt NUMBER r 12 * ?H * AREA= * F9. 3 t 8H t BURDEN) 

CALL HWS(XB) 

305 CONTINUE 
310 CONTINUE 

WRITE(6*6Q8) 

60B FORMAT (//SOH BURDEN TOTALS BY SURFACE ----------- -//) 

A ~ 0 o 

CALL HES ( 505 »K0 ) 

DO 320 J=l»4 
JJ=5U0+J 
CALL HAS(JJ) 

WRITE ( 6 1 606 ) JpAR(J) 

606 FORMAT (/eh SURFACE * 12 » 13H * TOTAL AREA-*F9.3r 14H* TOTAL BURDEN) 
CALL HWS(JJ) 

A=A+AR( j) 

320 CALL HCS(5Q5eJJ*505s 1) 

CALL HAS ( 505 ) 

WR ITE ( 6 * 609 ) 

609 FORMAT (//50H ------------------------ -//) 

WRITE ( 6 * 6Q7 ) A 

607 FORMAT C//28H TOTAL, t ALL SURFACES. AREA= * F9 . 3 * 9H > BURDEN) 

CALL HWS ( 505 ) 

JT(NT)=KT 

XMT(NT)=AVF(S05) 

XVT ( NT ) ~AVF ( LL ) 

IF (L4.LE.0) RETURN 

L4=L4+530 

CALL HES (L4 * 505 ) 

DR(L4f i)=A 

RETURN 

END 
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SUBROUTINE HCSMAr IB» IC»K) 

C 

C HISTOGRAM COMBINING SUBROUTINE 

DIMENSION CAU1) »CB(U) #CR(2U »ZR(2l) r QR < 11 * 11) » NQ < 11) 

COMMON KK » IK *KR # RUN ( 7 ) » KS # STG ( 7 ) t KT » TSK ( 7 ) »DSC(6) »XDs 

• NE » ALS ( 4 ) r ItC ( 10 ) » I tT C 10 ) » AEC (10) » AET ( 10 ) t AED ( 10 ) » AEF ( 10 »4) » 
♦NO# 101(20) # 100(20) HOC (20) # 

• NP » UAC ( 120 r 2 ) # 1U1120) » I AB ( 120 » 6 ) » AAG ( 120 ) » AAS ( 120 ) #AAT(120) # 
.ND * NX 1 572) #0R(572#11) rXR(572»ll) #K0# 

• N S » K N T (20) » IR ( 5 ) »AR(4) » APA ( 2 ) t 

.NT# JT (100 ) * XMT ( 100 ) »XVT(100) # Ll » L2 t L3 > L4 t KTS » iTE 
NA=NX ( IA) 

NQ=NX( IB) 

A-DR ( I A # 1 ) 

B=DR(IB# 1) 

X=XRUA#1) 

Y=XR(IB»1) 

CR ( 1 ) =2F l A # B » K ) 

CM=CR ( 1 ) 

QR(IC»1)=CRU) 

2R(l)=ZF(X»Y»K) 

IF ( NA . GT • 1 ) GO TO 50 
IF ( X • NE » 0 ♦ ) GO TO 30 
IF((K-3)*(K-4) »EG> « 0 ) GO TO 35 
30 IF ( N6 * GT « 1 ) GO TO 40 
35 XR(IC#1)=ZR(1) 

NX (10 = 1 
GO TO 1000 
40 NC=NB 

BO . 45 U=2 *NC 
Y=XR ( IB # J ) 

CR( J)=DR(-IB» J) 

45 ZR(J)=ZF(X»Y»K) ; - - & 

GO TO 60 

50 IF ( NB . GT « 1 ) GO TO 90 
NC=NA 

DO 55 J=2 1 NC 
X=XR ( I A » J ) 

CR(J)=DR( I A # J ) 

55 ZR( J)=ZF(X» Y#K) 

60 NX(IC)=NC 

IF(ZR(NC) .GT.ZR(l) ) 60 TO 70 
XR ( IC » 1 ) =ZR ( NC ) 

JC=NC 

DO 65 J=2tNC 
DR(IC» J)=CR(JC) 

JC=JC“1 

65 XR(IGi J)=ZR(JO 
GO TO 80 

70 XR(lOD=ZR(l) 

DO 75 J=2 1 NC 
DR ( IC »J ) =CR ( J ) 

75 XR( IC# J)=ZR( J) 

80 IF ( K • LE « 4 ) GO TO 1000 
C ELIMINATION OF DUPLICATE X VALUES FOR K=5 

CR(1)=0. 

C K ( 2 ) - 0 • 

JC = 1 

ZR ( 1 ) =XR ( IC # 1 ) 

DO 85 J— 2 1 NC 

ruin -ro l 1 1 4-riP r Tr ..i) 
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IF ( XR ( IC » J ) *LE • XR ( IC » J“1 ) ) GO TO 05 

UC-JC+l 

CRC JC)=CR( J) 

2R(JC)~XR( IC» J) 

85 CR( J+1)=CR( J) 

NC-JC 
GO TO 300 

90 IF ( K »EG * 5 ) GO TO 2ff0 
JA=0 
JO-O 
JC=0 
KA=NA+1 
KB-NB+1 

IF(NA.GE.NB>GO TO 94 

NC“NB 

KC=KB 

GO TO 100 

94 NCsNA 
KC=KA 

100 UC=JC+1 
KC-KC-l 

IF(JC.GT.KC>GO TO 140 
IF( JA+l.GE.KA)GO TO 124 
JA~ JA+ 1 
KA~KA“ , 1 
X1~XR ( I A > uA ) 

X2=XR( lAf-KA) 

IF(JB+l.GE.KB)GO TO 120 

UB—JB+l 

KB“K8~1 

YI=XR( IB» JB) 

Y2=XR( IB*KB) 

IF (K®£0a2)G0 TO HO 
21 -ZF (XIpYIpK) 

Z2=2F(X1#Y2#K> 

IF (Z2®GT ®Z1 5 60 TO 102 
ZR( JO-22 
GO TO 104 
102 2R( JC)=Z1 

104 Z1=ZF(X2» YlrK) 

Z2=ZF(X2e Y2eK) 

IF ( Z2 . GT . Zi ) GO TO 106 
ZR(KC)=Z1 
60 TO 100 
10b ZR ( KC 5 ~22 
GO TO 100 

110 ZR( JC)=2F(Xlr Y1»K> 

ZR(KC)=ZF (X2» Y2»K) 

IF(ZR(JC)-ZR(KC) )100, U2H12 
112 2R< JC)=.5*(ZR(JC>+2R(KC) ) 

114 KC-KC+i 

IF (KC «GT»NC ) GO TO 116 

JC= JC+1 

ZR ( JC ) ~Z R ( KC ) 

60 TO 114 
116 NC=JC 

GO TO 140 

120 Yl=#5*(XR(I8*vJB>+XRnB»KB) ) 

ZR< JC ) =ZP 1X1 * Y1 ?K ) 
ZR<KC)=ZF(X2»Y1»K> 

60 TO 100 



124 Xl=.S*(XR(IArJA)+XR(lA#KA) ) 

IF C JB+l.GE.KBJGO TO 128 
JB+ A 
KB=KB~1 
Y1=XRC IB» JB) 

Y2=XR( IB* KB) 

IF ( K «EQ ® 2 ) GO TO 150 
21=ZF (XI* Y1 *K) 

Z2.-ZF ( Xi e Y2 *> K ) 

IF ( 22 a GT « Z1 ) GO TO 126 
ZR( JC)=Z2 
ZR(KC)=Z1 
60 TO 10U 
126 ZR(JC)=2i 
ZR(KC)=Z2 
GO TO 100 

128 Yl=.5*(XR(lB*oB) + XR( IB * KB ) ) 

ZR( JC)=ZF(XirYl»K> 

ZR(KC)=2R( JO 
GO TO 100 

130 2R ( JC ) -2F IX1»Y2»K) 
ZR(KC)=ZF(X1*Y1*k) 

GO TO 100 
140 NX ( IC)=NC 

IF(NC.EQ.l)GO TO 35 
IF{K.EQ.2)G0 TO 171 
DO IbO J=1#NC 
150 CR(J)-0 9 

DO 170 JA=2»NA 
DO 170 JB-2 * NB 
P=DR(IA* JA)*OR(IB» JB) 

X1=XR(IA* JA-1) 

X2-XR ( I A » JA ) 

Y1=XR(IB»JB-1) 

Y2=XR< IB* JB) 

Z1=2F ( XI » Y1 p K ) 

22=2F(X1»Y2pK) 

IF ( Z2 . GT • 21 ) GO TO 154 

2A-Z2 

GO TO 155 

154 ZA-Zl 

155 21-ZF ( X2 e Y i e K ) 

Z2-ZF(X2?Y2*K) 

IF ( 22 9 GT o 21 ) GO TO 158 
ZB-Z1 
GO TO 160 
158 ZB=22 

160 00 170 JC=2»NC 

IF(ZR(JC) .LE *ZA ) 60 TO 170 
IF(2R( JO e Gt.2B)GO TO 166 
CR( JC)=CR( JC)+P*(ZR(uC)-Za)/(ZB-ZA) 
GO TO 170 

166 CRC JO=CR( JO+p 
170 CONTINUE 
60 TO 300 
171 KC~1 

DO 172 JC=2*NC 

1 F ( ZR ( JC ) . G T . ZR ( NC ) ) 2R ( JC ) =ZR ( NC ) 
XFCZR(KC) .GE.ZR(JC) ) GO TO 172 

KC-KC + 1 
ZR(KC)=ZR( JC) 
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172 CONTINUE 
NC-KC 

lF(NC.GTol)60 TO 17$ 

NX (10=1 
XR(IC»1)=2R(1) 

GO TO 1000 

173 DO 174 JA=2rNA 

1 f4 C A ( JA ) =DR ( X A » JA ) 

DO 100 JC=2»NC 
21=2R( JC-1) 

Z2.-ZR ( JC > 

DO 175 JA=1*NA 

175 GR< JA» JC)=0. 

DO 178 JB=2»NB 
Y1=XR(IB»JB-1) 

Y2=XR( ZBr j0> 

XA=2F(Zl?Yi»l)+.000i 
XA=2F ( Z1 » Y 1 » 1 ) 

X0~2F ( 22 f Y2 » 1 ) 

DO 178 JA~2 1 NA 

IF (XR< IA? JA) 6 LE.XA)G0 TO 178 

IF(XR(IA» JA) • GE • X8 ) GO TO 176 

QR ( JA » JC ) =QR ( JA » JC ) +DR ( IB » JB ) + ( XR ( I A » JA ) -XA ) / (XB-XA ) 
GO TO 178 

176 QR ( JA» JC ) =GR ( JA > JC ) +DR ( IB » JB ) 

178 CONTINUE 

KA=NA 

UO 180 JA=2 » NA 

OR(KAf JC)=QR(KA» JC)-QR(KA-1» JC) 

K A=KA~1 
180 CONTINUE 

DO 100 JA=2»NC 
KC~0 

QRX=o0001 
DO 182 JC=2»NC 

IF ( QR ( JA » JC ) .LE.0RX)G0 TO 182 
R C J C 
NQ(JA)=KC 
QRX=0R( JA* JC) 

182 CONTINUE 

DO 184 JG=2 ? NC 

184 QR( JAf JC)=QR( JA? JO/QRX 
CA(JA)=CA( JA)/GRX 
00 188 J=2 > NA 
IF ( J «EQ » JA ) GO TO 188 
QX=GR( J»KC) 

DO 186 JC=2#NC 

QR( J?JC)=QR( J? JC)-GX*QR( JA? JC) 

18b CONTINUE 

CA( J)=CA( J)”QX*CA( JA) 

188 CONTINUE 
190 CONTINUE 

DO 192 JA=2»NC 
KC=NQ( JA) 

192 CR(KC)=CA< JA) 

CR ( 1 ) = 0 e 

0 0 1 9 4 J A = 2 t N C 

1 F ( C R ( J A ) « L f & U <> ) 0 R ( J A ) ^ 0 $ 

194 CR(JA)-CR(JA)+CR(0A-1) 

GO TO 300 
200 JA=1 
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JB = 1 

JG=Q 

KASO 

KB-Q 

CA C 1 ) ■•* 0 » 

CB(X)=Q. 

NM-NA 

1 F ( NH « LT 9 NB ) NF'-NB 

205 IF(XK(IA» JA)-XR(IB» JB) ) 210 f 235 >255 
2X0 IF ( KA . GT s 0 ) GO TO 260 

IF(XR(IA»JA).LT.XR<IB»1))G0 TO 225 

DB-Cb ( JB“1 ) +(CB( JB>-CB( JB-1) ) * ( Xr ( I A » JA ) »XR ( IB t JB-1 ) ) / 
« (XR(ib> JB)-XR(If3i JB-1) ) 

GO TO 220 
215 UB=CB(NB) 

220 JC=JC+1 

CR( JC)-CA( JA)*DB 
2R( JC)=XK( IA» JA) 

225 IF ( JA.GE»NA) GO TO 250 
JA- JA+1 

CA( JA)=CA( JA-l).+OR(IA»JA) 

GO TO 205 * 

230 IF ( KB. GT * 0 ) GO TO 280 A 

KA=1 

GO TO- 205 
235 JC=JC4-1 

CR( JC)=CA( JA)*CB( JB) 

ZR( JC)=XR(IA» JA) 

1F( JA.GE.NAIGO TO 245 
JA=JA+1 

CA( JA)=CA( JA-1)+DR( IA»JA) 

240 IF ( JU.GE.NBJGO TO 250 
JB=JU+1 

C8( JB)=CB( JU-1)+DR(IB> JB) 

GO TO 205 

245 IF ( KB 8 GT e 0 ) GO TO 280 
KA=1 

GO TO 240 

250 IF ( KA . GT. 0 ) 60 TO 280 
KB=i 

GO TO 205 

255 IF ( KO « GT » 0 ) GO TO 215 

IF ( XR ( IB > JB ) #LT 8 XR ( IA» 1 ) ) 60 .TO 270 

DA=CA( JA-1)-MCA( JA)-CA( JA-l) )*(XR( IB> JB)-XR<1A» JA-1) )/ 
. (XR( IA> JA)-XR(IA# JA-1) ) 

60 TO 265 
260 DA-CA ( NA ) 

265 JC=JC+1 

CR( JC)=CB( JB)*DA 
ZR< JC)=XR( IBs JB) 

270 IF ( Jb.6E.NB)GC TO 275 
JB- JB+1 

CB( JB)=CB( JB“1)+DR( IB? JB) 

GO TO 205 

275 IF ( KA . GT 8 0 ) 60 TO 280 

Kli— 1 

GO TO 205 

280 IjC-JG , . 

290 IF (NC.LE.NM)GO TO 300 ’ ' 

0 J *"■ 1 ^ 

J J-2 
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NC=NC-1 
DO 294 J-2 » NC 

DS=(CK< J-1)+(CR( J+i)-CR( J-l) )*(ZR( J)-ZR(J-l) ) / <ZR ( J+l )-ZR ( J-l ) ) 
®-CR(J> >**2 



IF IDS. GE. DU) GO TO 
D J-DS 
J J“ J 

294 


IF(OJ.LT. .OOODSO 

TO 296 

294 

CONTINUE 


296 

DO 298 J=JJ»NC 
CR( J)=CR(J+1) 


298 

2R( J)=ZR( J+i) 
GO TO 290 


300 

F-l ® /CR ( NC ) 

00 302 J-2 f NC 


302 

CR< J)~F*CR(u) 
XR( ICrl)=ZR(l) 


308 

DO 310 J=2»NC 



tjRdCt J)=CR( J)-CR( J-l) 

310 

XR(IC»J)=2R( J) 
NX ( IC ) -NC 


1000 

DR(lCil)=AVF(IC) 
IF (K*EG«b) RETURN 



IF(DR(IC#1) .EQ.O. 
F=CM/DR(IC»1) 

CALL HMS(IC»FrIC) 

RETURN 

END 

5 RETURN 



SUBROUTINE HAS(N) • 

HISTOGRAM ADJUSTING SUBROUTINE 

COMMON KK» IK»KR»RUN(7) »KS»STG(7> ,KT.TSK(7) *DSCt6) . ID. 
.NE.A£5(4) ,IEC(10) .IET<10> .AEC(lO) »AET<10) .AED(IO) .AEF<10.4> . 
.NQ*I0T(2Q) »IOG<20) .I0C120) . 

.NP#DAC(120»2) »IUQ2Q) , I AB ( 120 . 6 ) . AA6 (120 > . AAS (120 ) » AAT ( 120 ) » 

. NO. NX ( 572 ) * DR ( 572 . 1 1 ) .XR(572»11> »K0» 

• NS»KKT(20) »IRC5) »A'R(4) »APA(2) » 

• NT. JTUQO) . XM 1 (1005 . XVT( 100) .Ll*L2»L3»L4.KTS» I'TE 
DIMENSION C(il> »D(U) .Z(ll) 

M=NX(N) 

IF(M9EQ.l) return 
IF(ID.LE.O) RETURN 
IF(NX( ID) .LE.O) RETURN 
AVER=AVF(N) 

MD=NX( 10) 

IF(MD.GT.l) GO TO 5 
XR(N»1)=AVER 
NX(N)=1 
RETURN 
5 C(1)=0. 

D(l )=09 

Z(MD)=XR(N»M) 

DO 10 I=2.M 
10 C(I)=C(I-1)+DR(N*I) 

IF (C(M) «EQ 9 1 • ) GO TO 19 
IF(C<M) .EQ.O.) RETURN 
DO 15 1=2 * M 
15 C ( I ) “C ( I ) /C ( M ) 

19 00 20 1=2 f MU 

20 D(I)=D(I“1)+DR(ID»I) 

IF ( 0 ( MD ) .EQ.l. ) GO TO 27 
IF(D<MD) .EQ.O. ) RETURN 
DO 25 1=2 . MD 
DR( ID. I )=DR( ID. I >/D(MD) 

25 D(I)=D(I)/D(MD) 

27 IF(D(MD-1) . N£ . 1 . ) GO TO 30 
MD=MD“1 
NX ( ID ) =MD 
GO TO 27 

30 1=2 
J=2 

31 QR(N*I>-OR(ID»I> 

32 IF(CU) .LE.D(I) ) GO TO 33 

2 ( I ) = XR ( N» U-l ) + (D ( I ) -C ( J-l ) )*<XR(N» u)-XR(N'J-l) ) / ( C C J) -C ( J-l ) ) 

1=14-1 

GO TO 31 

33 J=J+1 

IF( J.LE.M) GO TO 32 
DO 35 I=2»MD 
35 XR (N» I ) ~2 ( I ) 

NX ( N) =MD 
AVR=AVF (N) 

IF( AVR.fcQ.O. ) RETURN 
F=AVER/AVR 
DO 40 1=1. M 
40 XR(N»X>=XK(N*I>*F 
RETURN 
END 
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SUBROUTINE HES(J#K') 

HISTOGRAM EQUATING SUBROUTINE 
SET J=K, 

COMMON KK»XK*KR*RUN<7> *KS?STG(7) ,KT»TSK(7) rUSc( 6 ) * ID» 
e NE#AES(4) flECUO) rlET(lO) ?AEC(10) , AET(IO) » AED ( 10 > »AEF(10»4 ) t 
9 NO ? JOT C2G ) » ICQ (20 5 o ICC (20)* 

.NP»QACU2Q»2) #IU(12U) 0lAB(12G*6> eAA6(120) »AASC12Qj s»AATC 12Q) » 
.N0rNX(572) ? DR <572, 11) »XR<572, 11) »K0» 
o NS e KKT ( 20 ) , IR(5) # AR ( 4 ) , APA ( 2 ) » 

•NT»JT(1U0) »XKT(1U0) »XVT(100) , L 1 * L2 t L3 »L4 , KTS » JTE 
M=NX(K) 

NX( J)=M 

00 10 I- 1 , M 

DR(J»I>=OR<K»I> 

10 XR(J»I)=XR(K»I) 

RETURN 

END 
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SUBROUTINE HMS(IA#C# IC) 

r * 

v* 

C HISTOGRAM MULTIPLYING SUBROUTINE 

C HISTOGRAM IC EQUALS IA MULTIPLIED BY CONSTANT C 

COMMON KK# IK#KR»RUN(7) »KS#STG(7) »KT r TSK (7) #DSC(6) # ID» 

• NE » AES ( 4 ) » IcCdO) # IET(IO) »AEC(10) # AETdO) >AEOd0) #AEF(10#4) # 

.NO# I0T(20) » IOQ ( 20 ) » IOC (20) # 

•NP»DAC(120»2) » IU (120 ) »IAB(120»6> » AAG ( 120 ) » AAS ( 120 ) » AAT (120 ) # 

• ND » NX ( 572 ) #DR(572»U) #XR<572»11) »K0# 

•NS»KKT (20 ) » IR ( 5 ) » AR ( 4 ) »APA(2) # 

.NT# JT(IOO) » XMT ( 100 ) #XVT<100) #L1*L2#L3#L4#KTS» ITE 
IF (C.EQ.O. )G0 TO 20 
M=NX( IA) 

NX(IC)=M 
DO 10 J=1#M 
DR(IC# J)=0R(IA#U) 

10 XR(IC#J)=XR(IA# J)*C 
DR(IC#1)=DR(IA#1)*C 
RETURN 

20 NX(IC)=1 

DR(IC# 1)=0. 

XR(IC#1)=0. 

RETURN 

END 


SUBROUTINE HWSU) 

C HISTOGRAM WRITING ROUTINE 

C VERSION 16 JULY » 1969 & 37 

COMMON KK*XK»KR»RUN(7) »KSfSTG(7) »KT»TSK<7) »DSC(6) rXD» 

®NEe AES (4 ) , IEC ( 10) » XET ( 10) » AEC( 10) f ALT ( 10) r AEO( 10) eA£F< 10 »4) » 
s NOe IOT<20) eIQQ(20> # I0C(20) » 

»NP#DAC( 120*2) » IUU20) »IABU20»6) , AAG(12G) » AAS ( 120 ) t AAT C 120 > # 
•NO* NX ( 572 ) » OR (572 *11) » XR ( 572 * 1 I ) » KO * 

»NS t KKT ( 20 S ? IR C 5 ) *AR(4) » APA(2 ) , 

.NT* JTUOO) *XMT(iOO> * XVT ( 100 ) *Ll*L2*L3*L4*KTS»lT£ 

M=NX t I ) 

IF(M.LE.l) GO TO 100 
Q=AVF( I ) 

WRITE ( 6 * 600 ) U 

600 FORMAT (1H+*46X«13HMEAN VALUE r *E12.5) 

WRITE (6 *601 ) (OR( I* J) » J=2*M) 

601 FORMAT ( 15h PROBABILITY = * 5X * 10F 1 0, 5 ) 

WR XTE (6 * 602 ) ( XR < I » J ) * J=1 * M ) 

602 FORMAT ( 10H RANGE = »UE10.2) 

WRITE (6*609) 

609 FORMAT?/) 

RETURN 

100 WR XTE (6 * 603 ) XR(I»1) 

603 FORMAT (iH+*46X*2H= *El2.5/) 

RETURN 

end 


cno 


SUBROUTINE MAS ( I » T ) 


MICROBIAL ACCRETION SUBROUTINE ' 

I IS THE PART AFFECTED. T THE TIME OF UPDATE * 

COMMON KKe IK.KR.RUN(?> #KS.ST6(7) »KT.TSK(7) »DSC(6) » ID. 
«NE»ALS(4) »IEC(10)»IETd0> .AEC(IO) .AETCLO) »AED(10> »AEF(10»4> » 
•NO. JOT (20) . IOQ (20 ) » IOC (20 ) . 

• NP.OAC U2Q.2) . IUU2G) » 1 AB ( 120 » 6 ) * AAG (120 ) »AaS(120> »AAT(120) » 
«ND»NX ( 572) » DR (572. 11) »XR<572» I IS. KG. 

•NS . KKT (20 ) . IR(5) »AR<4) *APA<2) * 

•NT. JT(100> .XMTC10Q) #XVT(100) .Ll.L2.L3.L4.KTS. ITE 
OT=T»AAT ( D 
1F(DT.LE.0 6 ) return 
AAT(I)=T 
IE=IAB ( I » 5 ) 

IF( IE.LE.O ) RETURN 
IC=IEC(IE) 

Il=NO+l 

I2=W0+2 

I3=ND+3 

C=AAG(I)^AEC(IE)/DRUC.l) 

DO 40 J=l#4 

IF(IAWd.J) .EO.O)GO TO 40 
lF(AEF(It#J) .EQ.O.JGO TO 40 
IB=XAU(I»J> 

A=0RUB» 1) 

F=C*AEF( IE» J)*A 
CALL HMS(IC.F.II) 

DR(£lf 1)=AVF(11) 

V=AET( XE)*AES(U> 

F=EXP(“DT/V) 

CALL HMS(IB»F»I2) 

DR(I2.1)=AVF(I2) 

f=v*(i»-f> 

CALL HMS ( 1 1 . F . I 3 ) 

CALL hCS ( 12 . 13 . 16 . 1 ) 

DR(1B»1)=A 
40 CONTINUE 
RETURN 
END 


FUNCTION AVFU) 

AVERAGE VALUE FUNCTION 

COMMON KK* XK*KR*RUN(7) »KSpSTG(7) *kT*TSK(7) *uSC<6) »1D» 
b NE 9 AES ( 4 ) > IEC ( 10 ) »IET(10) »A£C(10) pAETQO) » AEQUO) » AEFC 10»4> 
. NO » 101(20) » 100(20) # IOC (20) * 

• NP * DAC ( 120 » ) r IU( 120) 1 1 AB ( 120 » 6 ) » A AG ( 120 ) pAAS(120) ?AATC120) 
•NO t NX ( 572 ) pURC572tll) pXR(572» 11) pKOp 
o NS r KKT (20 ) * IRIS) *ARU) . APA(2) » 

.NT#JT(100) pXMTUOO) »XVT llOO) *L1»L2>L3»L4#KTS»ITE 
M=NXl I ) 

IF(M.GT.l) GO TO 10 
AVF=XR ( I r 1 ) 

RETURN 
10 AVFsO • 

UO 20 

20 AVFsAVF+OR( I»0)*(XR( I» J-1)+XR(I» J) ) 

AVF=.5*AVF 

RETURN 

END 


non 


Function zf(x»y?k) 

2 FUNCTION OF X AND Y 

K determines the operation + - * /. 

GO TO Cx0?2Ge30#40e5G) ?K 
10 ZF=X*Y 
RETURN 
20 ZF-X-Y 
RETURN 
30 ZF-X^Y 
RETURN 
40 ZF"X/Y 
RETURN 
50 ZF=X 

IF ( Y » GT ® X ) ZF = Y 

RETURN 

END 


o r> o 


A4! 


SUBROUTINE UPS 

C DETAILED PRINTOUT SUBROUTINE (READS FROM TAPElD 

COMMON/X/M » MR * MG » MT , MS * DSC < 4 ) » N*NXi DR ( 11 > » XR< 11 ) »P< 5) »NP<5> 

REWIND 11 

READ (11) M6»MT» MSf (DSC ( J) » » 4 ) ? N » NX » (DR ( J ) » J”i t i 1 ) » - 

. (XR(J) »J=1» ID » IP(J> r J=l»5) » (NP(J) »J=1»5) 

WRITE ( b # 601 3 MR» C DR C j) rj=l»7) 

601 FORMAT (35H1MICR0BIAL BURDEN DETAILED PRINT0UT/4H RUN » 12 » 2X » ?A6// ) 
READ (5*500) K t NS r NT t NST » NS2 » NT2 * NST2 
500 FORMAT (715) 

k=o calls exit 

K-l PRINTS STAGE NS t TASK NT (SUBTASK NST) 

K-2 PRINTS FROM NS»NT»(NST) TO NS2»NT2» ( NST2 ) • 

10 READ ( 11 ) M»MR»M6»MT»MS» (DSC( J) » J=l*4) *N*NX» < DR < J) » J=l»li) t 
. (XR(U) * J=l»ll) t (P(J) *J=i*5) t (NP(J) * J=l»5) 

IF(M.LE.O) CALL EXIT 

11 IF(MG.NE.NS) GO TO 10 
IF(MT.NE.NT) GO TO 10 

IF< (NST.GT.O) .AND. (NST. NE. MS) ) SO TO 10 

IF (K»LQ. 2) GO TO 20 • 

15 CALL PLS 

READ (11) M * MR r MG * MT * MS r ( DSC ( J ) * J=1 # 4 ) » N * NX » ( DR ( J ) t J=l»ll ) t 
. (XR(U) » J=li 11) » (P( J) r J=i»5) » (NP(U) » J=l»5) 

IF(M.LE.O) CALL EXIT 
IF (MG .NE .NS ) GO TO 100 
IF(MT.NE.NT) GO TO 100 

IF ( (NST.GT.O) .AND. (NST.nE.MS) ) GO TO 100 
GO TO 15 
20 NS=NS2 
NT=NT2 
NST=NST2 
25 CALL PLS 

READ (ID M»MR*MG*MT »MS» (DSC ( J) #U=1 *4) *N*NXi> (DR< J) # J=l»ll) * 
a ( XR ( J ) t J~ loll) i» ( P ( J ) * J=1 o 5 ) 9 ( NP ( J ) » J- 1 » S ) 

IF(M.LE.O) CALL EXIT 
IF (MG.NE.NS) GO TO 25 
IF (MT.NE.NT) GO TO 25 

IF ( (NST.GT.O) .AND. (NST.NE.MS) ) GO TO 25 
100 READ (5*500) K t NS r NT » NST t NS2 » NT2 » NST2 
IF(K.LE.O) CALL EXIT 
GO TO 11 

END 



A42 


subroutine pl$ 

c parameter labeling subroutine 

C0MM0N/X/MeMR*MG»MT.MS.DSC(4) , N » NX t UR ( 1 1 ) , XR ( 1 1 ) » P C 5 ) , NP < 5 ) 

600 FORMAT(//40X»3'HRUN'»I2*7H» STAGE?I3?6H* TaSK»I4‘»9H» SUBTASK » 13 r 2X » 

« 4A6 ) 

GO T0(l*2»3»4»5»6»7?e*9) »M 

1 WRITE (6 #601 3 MR» (0R(U> » J=l»7 ) »MG‘» <DR< J) » J=8r U) » ( XR ( J) » J-1.3) 

601 FORMAT (35H1MICR08 1 AL BURDEN DETAILED PRINTOUT/4H RUN * 12 » 2X 1 1 kb/ 

®6H STAGE»l3r2X»7A6//.) 

RETURN 

2 WRITE (6 * 600 ) MR » MG »MT »MS r (DSC ( J>» J=1 r 4 ) 

WRITE (6.602 ) NP(1»NP(2> 

602 FORMAT (31H+ START OF NEW SUBTASK /23H PREREQUISITE SU8TASK 

•S“»I4»1H»I4»S2X. 10HSTART TIME) 

CALL HPS 
RETURN 

3 WRITE (6. 600 ) MR t MG >MT » MS t (DSC ( J ) » o=l r 4 ) 

WRITE (6. 603) NPU) »NP<2) »P(1) »NP(3) »P(3) #NP(4) .P(4) 

603 FORMAT (26H+PART BURDEN AFTER FALLOUT/5H PART. I4»7H» SURF. . 12, 

• 7H. AREA” . F7 . 2 r 7H . OPER . . 13 » IH » » F6 . 2 . 6H HOURS. 

• 10M» ENVIRON. . 13 . 1H. .F4 . 0 »5H MEN . » 10X » 6HBURDEN ) 

CALL HPS 

RETURN 

4 WRITE (6 » 600 ) MR t MG* MT . MS » (DSC ( J > . J=1 . 4 ) 

WRITE (6 t 604) NP ( 1 ) » NP ( 2 ) » ( P ( J ) » J=1 t 5 ) 

604 FORMAT (26H+PART BURDEN AFTER C0NTACT/5H PART » I4.7H. SURF. r 12. 

• 7H» AREA- . F7 » 2 » 1QH» CONTACTS . F5 . 3 » 5H . S1=»F5.3»5H* S2=f>F5.3* 

. 7H» T00L=.F7.1,7X.6hbURDEN) 

CALL HPS 
RETURN 

5 WRITE (6 .600 ) MR . MG » MT » MS . (DSC ( J ) » J=1 . 4 ) 

CALL HPS 

RETURN 

6 WRITE (6.600 ) MR t MG »MT t MS * (DSC ( J ) . J=1 » 4 ) 

RETURN 

1 WRITE C6.600 ) MR # MG # MT » MS » ( DSC ( J > t J=1 . 4 ) 

RETURN 

8 WRITE(6.G00 5 MR » MG t MT » MS » (DSC ( J ) » J=1 . 4 ) 

WRITE (6. 608) NP < 1 ) # ( P < J) » J=3» 4 ) * NP ( 2 ) # P ( 1 ) t P ( 2 ) »NP C 4 ) 

608 FORMAT (25H+PART BURDEN AFTER UP0ATE/5H PART. 14 . 2X » 2A6. ?H » SURF • t 
e 12 . 7H » AREA- . F7 . 2 . 7 h . TIME=.F8«3. 10H» EMV IRON . * 1 3 . 14X * 6HBURDEN > 

CALL HPS 
RETURN 

9 WRITE (6 » 600.) MR » MG # MT * MS . (DSC ( J > » J=1 . 4 ) 

WRITE (6. 609) NP (1 ) » NP ( 2 ) # P < 1 > . NP ( 3 > r P ( 4 ) . P ( 5 ) 

609 FORMAT (34H+PART BURDEN AFTER DECQNTAMINATIQN/5H PART » 14. 7H» SURF . » 

• I2»7H» AREA- . F7 . 2 . ?H . OPER . . 13 » lH . » F6. 2 » 25H HOURS. FRACTION REMOVE 

•0=.F6o3.8X. 6HBURDEN ) A 

CALL HPS 

RETURN 

END 




